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, 1.  INTRODUCTION 

i 


The  principles  upon  which  sound  ranging  is  based  were  formulated  in  the 
second  decade  of  the  century,  during  World  War  I and  remain  equally  the  basis 
for  the  techniques  as  they  exist  today.  The  fundamental  assumption  is  that  the  sound 
rays  proceed  spherically  outward  from  the  source  to  the  receivers,  moving  with  the 
local  speed  of  sound  as  modified  by  the  wind  velocity,  along  a shortest  pathler^th, 
regardless  of  terrain,  charocter  of  ground  surface,  and  stratification  of  the  atmosphere. 
The  field  procedures,  as  set  forth  in  FM  6-122  (1964),  are  designed  as  follows.  The 
measured  differences  in  arrival  times  are  recorded,  then  modified  successively  by  a 
temperature  correction,  a curvature  correction,  and  a wind  correction,  all  calculated 
from  graphical  devices. 

The  continued  use  of  sound  ranging  by  artillery  services  throughout  the  world 
for  more  than  a half-century  testify  to  the  essential  validity  of  these  assumptions 
and  to  the  techniques  adopted  for  their  field  applicotion.  However,  it  does  not  follow 
that  a reconsideration  of  the  physics  involved  cannot  lead  to  an  improvement  in  the 
accuracy  of  results  under  certain  conditions  of  atmosphere  and  terrain. 

It  was  therefore  the  philosophy  of  the  principal  investigator  to  attack  the 
entire  problem  of  sound  ranging  ob  initio — from  the  ground  up,  both  figuratively  arKi 
literally  (since  important  results  are  obtained  for  certain  types  of  ground  cover).  The 
following  questions  ore  thus  asked  and  answered  in  this  report: 
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(1)  How  valid  is  the  assumption  of  linear  propagation  of  a signal 
rising  from  the  firing  of  an  8-inch  howitzer? 

(2)  How  important  is  the  directional  effect  of  the  howitzer,  that  is, 
does  the  sound  propagate  as  a spherical  wave? 

(3)  To  construct  the  most  accurate  formulation  of  the  basic  sound 
ranging  formulas  and  program  them  for  digital  computer. 

(4)  What  is  the  effect  of  atmospheric  density  and  temperature  stratification? 

(5)  What  is  the  value  of  "ray  tracing"  in  sound  ranging? 

(6)  What  is  the  effect  of  an  absorbant  ground  surface,  such  os  ploughed 
earth  or  meadow-land? 

(7)  What  is  the  effect  of  terrain  obstacles,  in  particular,  hills  of  hundreds 
of  meters  in  elevation? 

(8)  Is  it  possible  to  use  any  of  the  signal  information  received,  other  than 
the  time  of  arrival  ? 

In, considering  and  answering  these  questions  we  have  used  digital  computation 
of  the  full  non-linear  equations  of  motion;  analytic  formulations  followed  by  digital 
evaluation  of  the  results,  and  data  from  the  Artillery  Meteorological  Comparisons  of 
the  Atmospheric  Sciences  Laboratory,  at  White  Sands  Missile  Range,  Oct/Nov,  1974. 

Throughout  the  work,  the  consideration  was  uppermost  to  incorporate  all 
results  in  operational  form  for  sound  ranging.  However,  this  is  not  always  possible 
at  the  present  time  because  of  the  lack  of  empirical  studies  of  sound  transmission  in  the 


3 


low  frequency  range  employed  by  current  sound  ranging  instrumentation. 

Readers  interested  in  results  rather  than  methods  may  turn  to  the  section 
titled  Conclusions  and  Recommendations. 


2.  A COMPLETE  AND  GEOMETRICALLY  ACCUR/JE 
FORMULATION  OF  THE  BASIC  SOUND-RANGING 
EQUATION 


The  sound  ranging  computational  routine  of  FM6-122  wos  devised  for  field 
use  before  the  development  of  small  electronic  programmable  calculators.  Thus  we 
shall  derive  all  formulas  in  this  report  under  the  assumption  that  such  calculators  will 
be  available  under  field  conditions.  The  calculating  routines  are  not  long  or  involved, 
and  the  necessary  programs  are  provided  in  this  report. 

First,  we  shall  derive  tv«>  alternative  formulations  of  the  basic  geometric 
formulas.  In  Figure  2. 1,  S is  the  sound-source,  R the  distance  to  the  microphone 
and  R+s^  and  R+Sg  the  distances  from  M^  and  M^  respectively.  The  three 
microphones  are  assumed  to  be  in  a straight  line,  separated  by  distances  d and  d . 

^ w 

Applying  the  law  of  cosines  to  the  traingles  SM  M and  SM  M , we  have,  respectively, 

I O ^ W 

= (R+Sg)^  + (d2+  - 2(R+S3)(d2+  d^)  cosij.^ 

2 2 2 

(R+Sj)  = *3)  ^dgCR+s^)  cos 


Elimination  of  cos  ij»  yields  the  formula  for  R, 


2R  = 


-''3‘3*S'‘2>'V‘'3> 


(2.2) 


If  we  now,  for  simplicity,  take  the  cose  in  which  the  microphones  are  equally  spaced. 


we  have  the  simpler  relation. 


2R  = 


= d,  = d / 

2 3 ' 


2 2 2 

2^2  - ^3 

*3-2*2 


(2.3o) 


TTiese  relations  for  R may  be  substituted  into  a low-of-cosines  (2. 1)  to  obtain  the 

corresponding  formulas  for  cos  «|)  Letting 

o 


'>’V'3/ 


we  have 


D cos  111  = 
3 


4 ^*3  ''  °’^3  ~ ^2^  V3  • 

OSj(2S3-  ij)  - djdjD  - 


and  for  the  case  of  equally  spaced  microphones. 


2d^(4s- 3,3)*  2.333(53 -Sj) 
2d  cos  »3  = 5 — 5 5 

‘'V3-^2-‘3-« 


(2.3b) 
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In  Hie  sequel  we  shall  assume  equally-spaced  microphones  so  that  the 
formulas  are  shorter.  There  is  no  principled  difficulty  of  deriving  any  of  the  formulas 
for  the  general  cose  of  arbitrary  spacing  or,  for  that  matter,  of  arbitrary,  rather  than 
straight  line,  configuration. 

A second  derivation  is  based  upon  the  geometric  theorem  that  the  locus  of 
a point,  the  difference  in  distance  of  which  from  two  given  points  is  constant,  is  a 
hyperbola,  with  the  two  points  as  foci,  thus  a single  pair  of  microphones  determine 
that  the  source  S lies  on  a certain  hyperbola,  and  a second  pair,  another  hyperbola. 
The  intersection  of  these  two  hyperbolas  uniquely  specifies  the  source  location. 

Consider  the  line  of  the  three  equally-spaced  microphones  to  be  the  x-axis, 
with  the  origin  located  midway  between  and  (Figure  2.2).  Thus  is 
located  at  x = ^2,  M at  x = di/2  and  M at  x = 34/2.  The  hyperbola  with  foci  at 

2 V 

( ± 4^2,  0)  and  with  constant  distance  s^  has  the  equation  for  and 


2 .2  2 

*2  - *2 
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The  corresponding  equation  for  M and  M is  obtained  by  applying  the  same  formula 

2 w 

at  the  point  x = -d,  that  is  by  letting  x'  = x +d  with 


..2 


d2-s-2 


4 


(Here  of  course  $'  = s.  - s in  terms  of  Figure  2,1  ).  We  have  two  equations  in 

w 2 
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two  unknowns,  and  by  elimination  we  obtain 


d ,..2  .2x 


(s-"  - s") 


I L 2 (.2  .2v  .2  /.2  .2v2 

$2  (d  - s-  ) 4s2  (d  -s'  ) 


(2.4a) 


where  the  minus  sign  is  used  if  is  closer  to  s than  (i.  e. , if  the  signal 
arrives  first  at  This  value  for  x is  then  substituted  into  the  formula 


(2.4b) 


to  obtain  the  coordinate  of  s normal  to  the  base  line. 

The  formulations  (2.3)  and  (2.4)  are  exactly  equivalent,  the  former 
giving  the  result  in  polar,  the  latter  in  cartesian,  coordinates.  We  shall  refer  to 
polar  coordinates  hereafter,  since  these  are  the  natural  expression  of  range  and 
azimuth. 

We  now  turn  to  the  problem  that  really  concerns  us  in  this  report:  how 
to  specify  the  distances  s and  s_  from  the  time-differences  of  the  signal  arrivals 
at  and  M^.  If  the  spherical  sound  wave  is  perfectly  reflected  by  the 

ground  and  is  not  diffracted,  scattered,  or  advected,  then 

s.  = C At.  i = 2,  3 

I I 

where  At.  is  the  time  difference  in  the  arrival  of  the  signal  at  M,  and  M. 

I *11-1 


(2.5) 
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To  the  extent  that  these  assumptions  are  invalid  in  any  particular  event,  or  that 
instrumental  error  effects  the  value  of  the  Af,,  the  basic  sound  ranging  formulas 
(2.3)  will  give  results  more  or  less  in  error. 

It  does  not  require  detailed  error  analysis  to  see  that  a small  percentual 
error  in  estimating  the  s.  will  lead  to  a small  percentual  error  in  calculating  cos  if 

I w 

by  (2.3b),  but  a large  percentual  error  in  calculating  R by  (2.3a).  To  take  a simple 
instance,  let 

d = 2km,  R = 7.07  km,  cos  »|i  =0.874 

3 

corresponding  to  true  values  of  s^  = 1.531  and  s^  = 3.225.  If  we  use  the  estimates 
s = 1,53  and  s = 3.25 — an  error  of  less  than  one  percent — we  obtain 

cosij>2  - 0.862  R -5.58 

an  error  of  about  one  percent  in  cos  if^  and  21  percent  in  R.  It  connot  improve  the 
accuracy  to  use  different  but  equivalent  formulas — for  example,  (2.4),  In  cartesian 
coordinates.  The  great  sensitivity  of  the  calculations  to  input  errors  is  inherent  in 
any  ranging  system  using  the  given  information  inputs,  to  wit,  the  signal -arrival  time 
differences. 

All  this  is  well  known  to  persons  acquainted  with  sound  ranging  or  any  other 
type  of  ranging.  It  is  restated  here  for  completness  and  emphasis.  There  ore  two 
approaches  to  reluction  of  error  that  will  be  considered  in  this  report: 

(1)  A more  consistent  and — it  is  hoped — occurate  formulation  of  the  effect  of  a 
wind  field;  and  (2)  a suggested  means  of  incomporating  the  effects  of  diffraction 
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by  each  of  a number  of  physical  sources:  atmospheric  temperature  gradients;  an 
absorbent  ground  surface;  and  surface  irregularities  such  os  hills.  In  each  case 
it  will  be  shown  that  under  certain  conditions  in  nature  significant  errors  could 
arise  by  foilure  to  take  account  of  the  accompanying  effects  in  a proper  way. 


Distortion  by  a constant  wind  field 


The  geometry  of  the  sound  ranging  problem  is  so  complicated  in  the  case  of 
even  a uniform  wind  that  we  must  begin  the  derivation  from  scratch  in  order  to  keep 
track  of  all  approximations.  Let  the  speed  of  the  wind  by  U,  directed  along  the 
x-axis,  without  loss  of  generality.  The  wave  front  is  then  given  by  the  equation 
(Whithorn,  1974,  pp.  254-259) 


(x  - Ut)^  + y^  = t^.  (2.6) 

We  may  rewrite  this  equation  in  terms  of  R,  the  distance  from  the  source: 

R^  - 2x  Ut  = (c^  - U^)  t^  (2.7) 

This  configuration  is  illustrated  in  Figure  2.3  for  the  two  points  and  ^2" 

The  observed  data  are  U,  9 , and  the  time  difference  *'2  “ •'p  Obviously  it  is 
necessary  to  make  simplifying  assumptions.  The  single  assumption  we  shall  make 
is  that  the  Mach  number  is  small,  that  is 


or 


2 2 
« c 


U V < R^. 


(2.8) 


g.1 
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No  approximations  aro  made  concerning  the  curvature.  Given  this  assumption,  we 
can  write  (2.7)  in  the  form 

c t 

* C.[l  ^Mi] 

If  we  now  apply  this  result  at  the  two  stations  and  M^,  we  have 

X X 

sj»Rj-R,  - r;'i] 


or  letting  and  using  the  first  order  approximation. 


we  have 


Ut  ± MR  , 


*2  " + Mix^-  x^) 


(2.10) 


$_  - cAt_  + Md  cos  0 

2 2 


(2.11) 


with  a similar  approximation  for  s^.  When  these  relations  are  substituted  into  the 
ranging  equations,  a result  is  obtained  valid  to  the  first  order  in  the  Mach  number 
with  no  approximatioru  as  to  the  curvature,  that  is,  no  assumptions  that  c(/R  or 
(x^  - Xj)/R  are  small.  It  would  be  possible  to  carry  out  the  approximations  to  the 


second  order  in  the  Mach  number,  or  to  work  implicit  equations  of  form  (2.7), 
However,  the  present  approximations  are  satisfoctory  in  olmost  all  cases  in  nature. 

In  any  event,  the  ranging  formulas  os  here  written  are  a consistent 
approximation.  They  are  programmed  in  FORTRAN,  with  a listing  provided  in 
Appendix  A.  It  should  be  emphasized  that  this  simple  program  is  well  within 
the  capability  of  current  hand-held  programmable  computers,  an  ideal  and 
flexible  piece  of  field  equipment. 
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3.  RAY  THEORY,  DIFFRACTION,  AND  SCAHERING 

U was  menHoned  in  Section  1 that  sound  ranging  is  based  upon  the  assumption 
that  a ray  travels  along  the  ground  surface  with  the  local  speed  of  sound.  It  is  the 
purpose  of  this  section  to  relate  this  assumption  to  the  theory  of  ray-tracir>g,  and  to 
the  theories  of  diffraction  and  scattering. 

The  fundamental  partial  differential  equations  for  families  of  sound  wave 
fronts  and  for  the  associated  ray  paths  were  set  forth  formally  by  Milne  (1921).  We 
mcy  gain  insight  into  their  structure  by  taking  a simple  instance  of  two-dimensional 
propagation,  say,  in  the  (x,z)-plane,  with  the  speed  of  sound  varying  only  in  the 
z-direction.  This  is,  of  course,  a reasonably  realistic  direction. 


The  ray-path  equations  are  then  (Milne,  1921,  p.  102) 


•j-  = tan  9, 
dx 


together  with  the  refraction  equation. 


c(z)  sin  9 = C sin  9 . 
o o 


(3.1a) 


(3.1b) 


Here  9 is  the  a.ngle  of  the  ray  with  the  vertical,  c(z)  is  the  stratified  sound  speed, 
and  the  subscript  zero  denotes  initial  values  for  the  ray  being  traced.  The  integral 
formulation  of  (3.1)  is 


X = X + sin  9 r 


c(z)dz 


z 2.2.2 
o c -sm  9 c 


This  is  consistent  with  our  qualitative  understanding  that  rays  are  deflected  concovely 
toward  lower  temperatures.  For  example,  if  c decreases  upward,  (3.2)  shows  that  for 


* 
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Figure  3. 1 (following  page)  Ray  Structure  for  Various  Atmospheric  Lapse  Rates 
of  Temperature  (after  Franck  and  Soger  (1974)). 

Abscissa  is  horizontal  direction  and  ordinate  is  vertical.  Plotted  lines 
are  rays,  labelled  in  degrees  of  initial  inclination  to  vertical.  Stippled  areas  ere 
zones  of  silence.  The  temperature  structures  in  each  of  the  cases  are: 

4a.  Lapse  rate  is  80  percent  dry  adiabatic. 

4b.  Temperature  inversion  surface  to  85  m. 

4c.  Temperature  inversion  surface  to  151  m. 

4d.  Temperature  inversion  surface  to  200  m. 
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0 > 0,  fhen  X and  z Increase  together,  that  is  the  ray  is  concave  upward.  Points  at 
o 

which  rays  are  oriented  horizontally  define  the  beginning  of  shadow  areas  into  which 
the  rays  given  by  (3. 2)  cannot  penetrate.  Such  ray  patterns  are  often  calculated  for 
atmospheric  conditions.  For  example,  Franck  end  Soger  (1975)  exhibit  their  ray 
calculations  (Figure  3. 1)  for  temperature  profiles  at  different  times  of  day,  indicating 
the  shadow  zones,  the  sound  source  being  at  ground  level. 

Calculations  such  as  these  are  correct  as  far  as  they  go;  but  they  must  not  be 
interpreted  as  presenting  a complete  picture  of  sound  transmission.  Formulations  of 
the  ray  equations,  whether  in  simplified  models  like  (3. 1)  or  in  terms  of  the  most 
general  geometric  optics,  describe  the  complete  problem  only  in  an  unlimited  medium. 
If  a boundary  is  present,  the  boundary  condition  becomes  part  of  the  problem,  and  the 
equations  must  be  reformulated.  The  shadow  boundary  is  defined  as  before,  but  now 
the  disturbance  propagates  over  the  boundary  into  the  shadow  with  the  local  sound 
speed.  The  wave  fronts  and  rays  in  the  shadow  area  are  said  to  be  diffracted.  This 
state  of  affairs  is  diagramatically  represented  in  Figure  3.2.  The  most  complete 
presentation  of  the  methods  of  geometrical  acoustics  is  that  of  Friedlander  (1958). 

Thus  the  basic  sound  ranging  assumption  is  in  a strict  sense  more  valid  than 
the  limited  ray  theory  such  as  that  developed  in  Milne  (1921)  and  applied  by  Franck 
and  Sager  (1975).  However,  despite  its  correctness  as  concerns  arrival  time,  this 
assumption  gives  no  Indication  of  the  amplitude  signature  of  the  arriving  signal.  W e 
shall  see  that  diffraction  drastically  alters  the  wave  or  pulse  form,  and  that  this 
alteration  may  produce  significant  errors  if  ignored  in  sound  ranging  procedures. 
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Whaf’  has  been  said  above  of  a stratified  atmosphere  is  equally  applicable 
to  the  problem  of  propagation  over  an  obstacle  or  over  an  absorbing  surface.  A trans- 
mitted pulse  propagates  into  the  geometric  shadow  as  a diffracted  pulse.  Although 
the  physical  mechanisms  are  similar,  there  are  quantitative  differences,  and  calculations 
must  be  performed  separately  for  each  physical  model. 

There  is,  of  course,  another  fundamental  objection  to  simple  ray-tracing,  which 
is  applicc^le  also  to  geometric  acoustics  in  general.  If  the  scotterer — whether  it  be  a 
temperature  gradient  or  a solid  obstacle — has  a characteristic  length  of  the  order  of  the 
wave  length  of  the  scattered  signal,  geometric  acoustics  becomes  a poor  or  invalid 
approximation  to  the  true  solution.  It  is  possible  to  treat  this  approximation  as  a first 
one  and  then  to  modify  it.  However,  for  purposes  of  this  study,  we  shall  drop  entirely 
the  physical  model  of  rays  and  ray-tracing  and  consider  each  cose  as  a boundary-value 
problem  of  mathematical  physics. 
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4.  TVIE  EFFECT  OF  AN  ABSORBANT  GROUND  SURFACE,  OR 
PLANE  vs.  SPHERICAL  WAVES  AND  PULSES 

When  a sound  source  radiates  in  the  atmosphere,  the  sound  travels  outward 
radially  as  a spherical  wave  or  pulse.  After  this  signal  has  traveled  for  some  minimum 
time,  the  distance  from  the  source  will  be  large  with  respect  to  any  length  scale  involved 
in  the  problem.  At  this  point,  one  might  be  tempted  to  ask,  can  the  model  adopted  be 
that  of  a plane  wave? 

A first  response  could  be  that  spherical  and  plane  waves  have  very  different 
dynamics,  one  amplitude  falling  off  inversely  with  the  distance  and  the  other  maintaining 
its  amplitude  unchanged  as  it  propagates. 

However,  this  difference  of  behavior  is  presently  of  little  importance  in  sound 
ranging,  and  we  turn  to  consider  a much  more  fundamental  objection  to  the  plane  wave 
model.  In  any  problem  of  artillery  sound  ranging,  the  presence  of  the  boundary  surface 
of  the  earth  is  a fundamental  consideration,  and  we  must  be  careful  to  treat  properly  the 
problem  of  propagation  near  a solid  boundary.  It  is  an  elementary  result  of  acoustics 
(Morse  and  Ingard,  1968,  p.  262)  that  the  reflection  coefficient  of  a plane  wave  at  a 
solid  boundary  (called  "locally  reacting") 


z cos9-  1 


(4.1) 


r z cos9  + 1 

where  z is  the  relative  impedance  of  the  boundary,  pc  that  of  the  air,  and  9 the  zenith 


angle  of  incidence.  If  z is  finite  and  ® ~ 2 ' ~ reflected  wove 

exactly  cancels  the  incident  wove.  A plane  wave  cannot  propagate  along  a boundary 
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of  finite  impedance.  On  the  other  hand,  if  z is  infinite,  then  c = 1 for  all  angles 
of  incidence,  and  the  reflected  wave  is  in  phase  with  the  incident  wove,  doubling 
the  amplitude. 

The  only  way  out  of  this  dilemma  is  to  adopt  an  adequate  physical  model: 
that  of  a spherical  wave  proceeding  from  a point  source  and  interacting  with  the 
boundary  of  arbitrary  impedance.  The  spherical  wave  is  mathematically  decomposable 
into  a continuous  spectrum  of  plane  waves  impinging  upon  the  plane  with  every  angle  of 
incidence  0 . Each  of  these  plane  waves  is  reflected  with  the  coefficient  given  by  (4. 1); 
so  that  when  these  are  summed  over  all  0,  the  resultant  wave  is  capable  of  propogation 
alor^  a plane  boundary  without  cancelation.  When  the  impedance  x is  very  large  in 
magnitude,  the  sphericol  wove  propagates  almost  undistorted,  except  that  its  amplitude 
is  doubled  by  reflection. 

It  should  be  mentioned  that  another  type  of  solid  boundary — in  addition  to  the 
locally  reacting — is  physically  possible,  called  a boundary  of  “extended  reaction". 

This  is  analogous  to  a two-media  system  in  which  the  ground  is  interpeted  as  a medium 
that  responds  as  an  acoustic  system  to  the  passage  of  the  atmospheric  wave.  The 
coefficient  of  reflection  in  this  case  is  given  by  (Morse  and  Ingard,  1968,  pp.  266  ff) 

m cos9  - 

c = 

m COS0  + 

where  m = p^/p  and  n = c/Cy  the  subscripted  quantities  referring  to  the  bounding 
medium.  The  locally  reacting  surface  leading  to  (4. 1),  we  interpret  as  one  which 
permits  penetration  by  air  motions,  but  does  not  transmit  horizontally  an  acoustic  wave. 
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This  may  be  realized,  for  example,  by  a ploughed  field,  or  a meadow  in  which  the 
I soil  is  kept  loose  by  the  grassy  growth. 

i The  treatment  of  the  problem  of  a spherical  wave  propagating  over  a plane 

i 

I surface  of  finite  impedance  has  a long  history.  It  was  first  solved  by  Rudnick  (1947), 

f \vho  applied  the  Sommerfeld  solution  for  electromagnetic  waves  to  the  acoustic  case 

I of  a surface  of  extended  reaction.  Loter  Lawhead  and  Rudnick  (1951)  and  Ingord 

(1951)  investigated  also  the  same  problem  for  a locally-reacting  surface.  These  and 
other  early  studies  are  reported  in  textbook  form  by  Brekhovskikh  (1960,  Ch.  4). 

Interest  in  the  subject  has  recently  revived,  and  a new  burst  of  research 
activity  has  introduced  new  methods  and  increased  our  detailed  knowledge  of  the 
solutions — see,  for  example,  Delaney  and  Bazley  (1S70);  Wenzel  (1974);  Chien  and 
Soroka  (1975);  Thomasson  (1976);  and  Donato  (1976a,  1976b). 

The  points  of  emphasis  here  are  (1)  that  especially  when  dealing  with  regions 
near  a boundary,  the  plane  wave  approxir»ation  may  be  inadequate  and  misleading, 
and  (2)  that  models  involving  spherical  wave  propagation  near  a boundary  quickly 
lead  to  problems  at  or  beyond  the  current  state  of  the  art  in  mathematical  analysis. 

In  the  ensuing  treatment,  we  shall  take  core  that  our  approximations  are  valid  physically 
as  well  as  mathematically. 

The  analyses  referred  to  above  have  in  common  that  they  apply  to  monochromatic 
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radiation  from  an  acoustic  source.  For  our  sound  ranging  problem,  a much  more  realistic 
model  is  the  pulse,  involving  the  complete  spectrum  of  waves.  This  is  porticulorly 
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important  given  the  fundamental  role  of  the  signal  arrival  time  in  current  techniques. 
We  shall  thus  restrict  our  attention  to  pulse  propagation. 

The  basic  analytical  work  here  was  performed  by  Doak  (1952).  It  differs  from 
the  techniques  in  the  monochromatic  case  by  treating  an  initial-value  problem  and 
studying  the  signature  of  the  signal  at  various  distances  from  the  source  as  a function 
of  time. 

We  assume  a simple  unit  step-function  pressure  pulse  in  three  dimensions 

_ H(t  - R/c) 

P - -JJR 

where  H(t)  is  the  heaviside  step-function, 

H(t)  = 0 t<0 

(4.3) 

= 1 , t>0 

and^  is  the  distance  from  the  source.  The  source  is  assumed  to  be  an  arbitrary  distance 
above  a uniform  locally-reacting  horizontal  plane  of  finite  impedance  z,  so  that  the 
reflection  coefficient  c for  each  component  wave  in  the  pulse  (4.2)  is  given  by 
(4. 1).  The  mathematics  involves  a number  of  coordinate  transformations,  but  is 
entirely  straightforward.  For  details,  the  reader  is  referred  to  Doak's  (195?)  paper. 

The  general  result  may  be  written  in  the  form 


H(t-  IVc)  ^ H(t-  r/c) 
411 R 4Tr  r 


’ ^ 

where  t = ct/r  and  F(9,  X)  is  the  inverse  Laplace  transform  representing  the  pulse  as 
diffracted  by  the  absorbing  boundary: 
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F(9,X)  = ^ J exp  if- ’^f  ^ ds  (4.5) 

[(l+zXcosG)  +2  (X  (X  - l)sin  9j 

Here  B is  the  Bromwich  path  in  the  complex  plane,  with  all  singularities  of  the  integrand 
on  i-he  left. 

The  geometry  of  (4.4)  is  given  in  Figure  4.1,  in  full  generality.  For  current 
sound  ranging  techniques,  S coincides  with  S'  at  the  surface,  and  M is  also  at  the  surface, 
so  that  R = r and  9 = ir/2. 

The  evaluation  of  (4.5)  depends  on  knowing  the  impedance  z as  a function  of 
wave-frequency,  that  is,  in  this  formulation,  as  a function  of  s.  As  a first  approximation, 
we  may  make  an  evaluation  for  z = constant.  In  this  case. 


F(9,X)  = 


6 (t  - rX/c) 

[(1  -hz  Xcos9  )^  -hz^(X^-  1)  sin^9  ) 


where  6 (t)  is  the  Dirac  delta-function. 

Substituting  this  expression  into  (4.4),  we  obtain 


H(t-R/c)  ^ H(t-r/c) 
4ir  R 4w  r 


[(1  +ZT  cos9)^  + z^(t^-  1)  sin^9  ]"  ^ 

4Tr  r 


where,  as  before,  = ct/1?. 
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The  first  two  terms  ore  the  incident  and  the  reflected  pulses, 
respectively.  The  third  term  shows  the  effect  of  the  locally-reacting  (absorbant 
surface.  The  arrival  times  ore  those  given  by  geometrical  acoustics,  but  the  amplitude 
is  distorted  by  the  diffracting  effect  of  the  surface  on  the  reflected  pulse.  Just  after 
arrival  (t  = 1),  we  have 


p = — ~ + — 1—  I 1 i i 

4ttR  4Trr  l+zcos9  ^ ' 

so  that  for  low  impedance  and  distant  points  ( G-*it/2),  the  diffraction  can  more  than 
cancel  the  reflected  pulse  amplitude.  When  the  signal  arrival  is  long  past  (t  » J), 
we  have 


that  is,  the  pressure  is  approximately  the  sum  of  the  incident  and  reflected  pulses. 

For  a more  detailed,  quantitative  analysis,  we  consider  the  special  sound 
ranging  case  for  which  0 = 11/2,  R = r.  Thus 


p = (4.7) 

4it  R 4it  R 

Incident  and  reflected  pulses  arrive  simultaneously  at  t = 1,  but  at  this  moment  the 
diffracted  pulse  exactly  concels  the  sum  of  the  incident  and  reflected  pulses,  so  that 
the  amplitude  is  zero.  The  diffraction  dies  out  as  time  increases,  but  there  is  a 
non-zero  time  interval  between  the  arrival  time  and  the  time  that  the  microphone  detects 
the  signol.  The  magnitude  of  this  interval  will  depend  on  the  impedance  of  the  ground 
surface  z on  the  distance  R of  the  microphone  from  the  source,  and,  of  course,  on  the 
sensitiveity  of  the  microphone. 
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Figure  4.2  Signature  of  Spherical  Pressure  Pulse  on  a Surface  of  Finite  Impedance 


The  cbscissa  shows  scaled  time  after  arrival  (ct/R  - 1)  and  the  ordinate, 
scaled  pressure.  The  curves  ore  labeled  with  the  relative  impedance  values  Z. 
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In  all  sound-ranging  theory  to-date,  it  has  been  assumed  that  the  signal 
arrival  is  a clearly  defined  event,  idealized  as  a step  function,  and  not  a gradual 
increase  from  a zero  amplitude. 

The  signature  (4.7)  has  been  plotted  for  a range  of  impedance  values  in 
Figure  4.  2.  It  v/ill  be  noted  that  the  behavior  is  strongly  dependent  on  the  impedance. 
For  z = 100,  the  curve  is  close  to  a step-function,  and  the  signal  rises  rapidly  to  its 
maximum  steady  value  of  two.  For  z of  the  order  of  unity  its  rise  is  very  much  slower 
and  the  signature  bears  no  resemblance  to  a step  function.  It  follows  that  a much 
longer  time  interval  (t  - 1)  would  pass  before  the  signal  rose  to  detectable  level. 

We  may  apply  these  results  to  the  sound  ranging  problem  as  follows.  Consider 
that  a certain  critical  pressure  p^  ;$  known,  which  will  activate  the  microphone  system. 
This  will  be  the  same  for  and  for  M^.  Thus  from  (4.7),  we  obtain 


% - r 2 2,2  p2  / 2 

X-R,  -[zct,-R,  fe-in 

_ p-1  , 2 2 2 p2  , 2 

= R2  ■ L2  c fj  “ ^2  * " ^ 


(4.8) 


where  X is  a known  constant.  Here  t^  and  t^  are  not  the  true  arrival  times,  but 
rather  the  times  at  which  the  signal  first  reaches  a detectable  amplitude  at  and 
respectively.  What  we  observe  therefore  is,  as  before,  the  time  difference 
tj  “ t^,  and  from  it  we  wish  to  deduce  the  distance  difference  R^-  Rj. 

A simple  mathematical  case,  although  extreme  as  a physical  assumption,  is 
that  for  z = 1.  Then  (4.8)  becomes 


R2“R,  =C(t2-f^)[(l-XR^)0  nctj)"']  (4.9) 


Fig.  4,3  Impeckince  of  a Grassy  Surface  for  Various  Frequencies  (after  Embleton,  Piercy,  and  Olson  (1976)). 
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In  (4.9),  the  bracketed  factor  on  the  right  represents  the  correction  due  to  the  low 
impedance.  The  classical  formula  AR  = cAt  is  recovered  by  taking  X -•  0,  that  is, 
infinitely  sensitive  microphones. 

Field  Studies 

Unfortunately  there  are  no  direct  field  studies  known  to  the  author  of  sound 
ranging  over  a surface  of  low  impedance.  However,  problems  of  aircraft  noise  have 
led  to  some  recent  Investigations  of  interest  on  acoustic  wave  under  realistic  outdoor 
conditions,  for  example,  Dickenson  and  Dock  (1970);  Piercy  and  Embleton  (1974); 
Embleton,  Plorcy,  and  Olson  (1976).  These  field  studies  concentrated  on  noise  in 
the  relatively  high  frequency  range,  of  the  order  of  10  and  10  Hertz.  Some  of 
the  results  are  reproduced  in  Figures  4.3  and  4.4.  The  curves  show  a tendency  for 
the  impedance  to  increase  with  decreasing  vwive  frequency,  but  there  are  not  suf- 
ficient experiments  at  the  lower  frequencies  to  permit  us  to  extrapolate  the  curves. 
Further,  the  low  impedances  were  observed  over  meadowland.  It  is  possible  that  a 
looser  surface — plowed  fields,  forest  ground  cover — could  reduce  the  impedance  still 
further. 

Propogation  of  an  N-wove 


The  striking  results  from  these  experiments  suggested  a more  realistic 
analytical  model.  The  solution  (4./)  is  for  a simple  step-function  pulse;  this  is  not, 
however,  the  typical  signal  form  projected  by  an  artillery  blast.  This  question  will 


be  discussed  in  detail  below,  but  we  can  here  use  the  simple  model  of  an  N-wove 
for  the  early  signal  form: 

N(t)  = 0 0 < t 


= l-y  0<t<T  (4.10) 

= 0 T<t 


The  solution  (4.6)  can  now  be  interpreted  as  an  influence  function  for  the  more  general 
solution, 

4irp  = 2N(t-R/c)-2  r*  N‘(t,)  H(t-^ -t.) 

'o  J dt  , (4.11) 

{l+z2[^(t-tj2+lf 

where  N'(t)  is  the  derivative  of  (4. 10),  and  we  have  assumed  propagation  along  the 
ground  surface  ( 0 = ir/2).  If  the  integral  in  (4.11)  is  evaluated,  we  have 


4TrRp  = 2N(t-R/c)  - 


2H(t-R/c) 


2H(t-IVc-p 


[1  +z2(t2-1)  [1  + z2(t-tV- 1]^ 


2 2 

4 , T T - a 

— ^ 
zT  1 + z 


,1  +22 

4 „ T + T - a 

★ ^ if  if  0 9 

zT  t-T+  (t-T  ) -a'' 


1 <T  < 1 + T 


★ 

1 + T < T 


(4.12) 


if  5-2 

Here  t = ct/R,  as  before,  T = cT/R,  and  a = 1-z  . In  (4. 12)  the  first  term 


represents,  of  course,  the  undistorted  pulse  traveling  with  speed  c,  the  second  and 
third  terms  represent  the  distortion  by  the  absorbant  boundary  of  the  impulses  at  t - 0 
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and  t = T,  respectively,  and  the  fourth  term  represents  the  distortion  of  the  linear 
fall-off  of  the  signal  during  the  period  0 < t <T. 

One  would  anticipate  that  surface  absorption  would  hove  a more  pronounced 
effect  on  a signal  of  short  duration  than  on  the  single  step  function,  and  this  is  the 
case.  Some  calculations  from  (4. 12)  are  graphed  in  Figure  4.5.  Here  we  see  that 
the  z = I curve  is  scarcely  distinguishable  from  the  axis:  the  N-wave  signal  has  been 
effectively  wiped  out.  The  curve  for  z = 100,  on  the  other  hand,  resembles  the  initial 
disturbance  very  closely.  This  impedance  is  not  appreciably  different  from  total 
reflection  in  form,  although  the  amplitude  is  reduced  by  about  30  per  cent. 

The  formula  (4.  12)  is  much  more  complicated  than  (4.7),  and  implicit 
formulation  for  the  sound  ranging  problem  corresponding  to  (4.y)  would  not  contribute 
to  our  intuitive  understanding,  and  hence  is  not  written  down  here.  We  appear  to 
be,  at  this  point,  very  far  from  the  original  sound  ranging  formulation;  and  this 
situation  may  be  described  as  follows. 

So  long  as  the  pulse  is  propagating  undiffracted,  it  must  arrive  in  advance 
of  any  reflected  signals,  and  if  it  is  of  detectable  amplitude  at  all,  it  will  register 
at  the  microphone  as  a pulse,  presumably  in  the  form  of  an  N-wave.  However,  once 
the  pulse  travels  as  a diffracted  signal  in  the  shadow  zone,  its  detectable  arrival  times 
no  longer  satisfy  the  fundamental  assumptions  of  sound  ranging.  The  problem  of  sound 
ranging  then  becomes  essentially  an  "inverse  problem"  of  deducing  local  information 
at  remote  sites  from  information  received  at  three  or  more  stations.  This  formulation 
will  be  discussed  in  more  detail  in  a later  section  of  this  report.  The  emphasis  of  this 
section  is  that  such  a problem  can  arise  even  if  the  medium  is  homogeneous  and  the 
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ferrain  is  level:  a type  of  diffraction  occurs  v/henever  the  earth  surface  has  a 
relative  impedance  of  the  order  of,  say,  ten  or  less. 


In  problems  of  diffraction  and  scattering  that  we  shall  consider  here  and  below, 
includir^  that  of  a stratified  atmosphere,  have  this  in  common:  they  are  too  difficult 
analytically  for  a complete  solution.  However,  fortunately  our  interest  is  in  the  time 
interval  soon  after  arrival  time,  and  for  this  brief  period  an  asymptotic  result  is  possible. 

This  result  is  in  a form  amenable  to  digital  evaluation,  and  such  computations  have  been 
performed  and  are  represented  in  graphical  form. 

As  was  pointed  out  in  Section  3,  simple  ray  theory  is  not  applicable  to  a 
stratified  atmosphere  unless  diffracted  rays,  propagating  into  the  shadow  zone,  are  token 
into  account,  and  the  amplitude  is  calculated  along  the  ray. 

The  proper  diogramatical  representation  of  the  ray  pattern  is  that  of  Figure 
This  is  drawn,  obviously,  for  an  atmosphere  In  which  the  temperature,  and  so  the  speed 
of  sound,  decreases  with  height.  Thus  the  rays  are  concave  upward.  That  particular  ray 
that  is  tangent  to  the  bounding  surface  at  its  lowest  point  defines  the  shadow  boundary. 

This  figure  may  now  be  compared  to  Figure  3.  la.  In  3. 1,X  fhe  source  is  on  the  boundary, 
and  therefore  the  entire  surface  z = 0 is  in  shadow.  According  to  scattering  theory,  the 
diffracted  ray  proceeds  at  sound  speed  into  the  shadow  zone,  which  it  is  thus  misleading 
to  refer  to  as  a "silent  zone".  The  problem  now  becomes  a classical  eigenfunction 
expansion. 

For  a simple  atmospheric  model  in  which  temperature  decreases  with  height,  a 
pulse-propagation  solution  was  worked  out  some  years  ago  by  Friedlonder  (1958).  It  is 
assumed  that  the  bounding  surface  is  perfectly  reflective. 


The  Green’s  function  for  the  problem  is  the  solution  of  the  equation  for  a 
step  function  impulse  at  the  origin 

2 

> . ^2  G = 6(x)  6(y)  6(z)  H(t) 

c 9t^ 

where  t>  is  the  delto  function,  with  the  boundary  condition  9G/92  = 0 at  z = 0. 
Friedlander’s  technique  is  to  take  the  Laploce  transform  of  this  equation  and  determine 
its  behavior  as  the  transform  variable  s This  gives  the  form  of  the  signal  for  small 

times  after  arrival,  a very  suitable  asymptotic  solution  for  our  purposes.  The  analysis 
is  straightforward,  leading  to  the  result,  for  z 0, 

4TrRp  = H(t  - 1)  exp  [ - ^ J,  (5.1) 

o 

where  t = c t/R,  as  before,  and  F = - (dT/dz)  is  the  lapse  rate  at  the  surface.  The 
o o 

form  of  (5. 1)  is  a special  case  of  a more  general  form  of  diffracted  time-signature, 

Rp  = A(t-1)"^  exp  [ - ] . (5.2) 

The  signal  starts  at  the  arrival  time  t = 1 with  zero  amplitude  and  grows  slowly  or 
rapidly  according  to  the  coefficient  B in  the  exponential.  In  the  more  general  case, 
with  the  exponent  q different  from  zero,  there  is  an  effect  tending  to  decrease  the 
amplitude  with  increasing  time  t-1.  However,  it  must  be  remembered  that  expressions 
like  (5.  I)  and  (5.2)  are  derived  as  asymptotic  expansions  for  (t-J)  small,  and  thus  must 
not  be  expected  to  give  reliable  results  for  long  times. 
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Some  signafures  have  been  calculated  from  (5.  1)  and  are  presented  in 
Figure  5.2.  We  have  taken  the  distance  from  the  source  as  3.3  km  and  plotted  two 
curves,  one  for  a lapse  rate  of  10  deg.'^km  and  one  for  5 deg/km. 

It  will  be  noted  that  both  signatures — unlike  those  of  Figures  y.  2 and 
y.  S for  low  surface  impedance — rise  rapidly,  resembling  closely  the  original  impulse 
H(t  - iVc).  It  is  evident  from  (5.  l)  that  if  we  were  to  choose  a very  large  temperature 
gradient  F,  or  large  distance  R,  it  would  be  possible  to  obtain  signatures  relatively 
flat  at  the  origin.  However,  within  the  practical  limits  of  sound  ranging  in  the 
atmosphere,  we  may  expect  Figure  5.2  to  be  typical. 

We  thus  conclude  that  it  is  highly  unlikely  that  normal  atmospheric  temperature 
gradients  under  normal  sound  ranging  conditions  will  contribute  significantly  to  the  ranging 
error.  This  conclusion  is,  of  course,  subject  to  revision  by  a study  of  instrumental  sensitivity. 
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6.  THE  EFFECT  OF  A TERRAIN  OBSTACLE 

l^  is  evident  that  irregularities  in  the  surface  over  which  the  sound  is  transmitted 
will  hove  an  important  effect  on  the  signal,  and  considerable  effort  has  been  devoted  to 
this  topic.  Reflections  from  periodically  corrugated  surfaces  are  discussed  by  Meechom 
(1956)  and  by  Proud  et  al.  (1958);  and  a more  general  method  of  treatment  is  suggested 
by  Meecham  (1956), 

However,  we  anticipate  that  the  most  satisfactory  method  of  treating  surface 
irregularities  will  be  statistical,  that  is,  by  statistically  parameterizing  the  properties 
of  the  surface.  This  can  be  done  relatively  simply  for  a single  plane  wave  when  the 
surface  has  an  impedance  or  roughness  witfi  known  variance  (Morse  and  Ingord,  1968, 
pp.  441-449).  We  are  not  aware  of  a corresponding  treatment  for  a spherical  ‘pulse. 

Thusthe  content  of  this  section  will  be  limited  to  terrain  obstacles  of  the  order 
of  the  wave  length  of  the  sound  radiation,  or  greater.  For  the  sound-ranging  problem, 
the  obstacle  considered  will  be  of  the  size  of  a hill  or  small  mountain.  Although  the 
mathematics  is  difficult,  the  analysis  can  be  done  by  the  approximations  and  techniques 
outlined  in  Section  5. 

Physically,  the  geometric  optics  of  this  problem  is  similar  to  that  of  a stratified 
atmosphere,  in  that  our  greatest  interest  lies  in  the  shadow  region  in  which  only  diffracted 
rays  can  penetrate.  The  picture  is  that  of  Figure  6. 1.  Again  we  shall  find  that  the  signal 
arrives  with  the  local  speed  of  sound,  but  with  zero  amplitude,  so  that  time  of  arrival 
and  time  of  detection  are  not  identical. 
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For  reasons  explained  more  fully  in  preceding  sections,  the  existence  of  the  earth 
as  a boundary  along  which  the  pulse  must  propagate  precludes  the  use  of  any  plane  wave 
model.  Thus  application  of  the  great  body  of  analytical  results  on  the  scattering  of  plane 
waves  by  obstacles  of  various  configurations  is  denied  to  us.  The  simplest  relevant  formal 
problem  is  that  of  a spherical  pulse  scattered  by  a cyclinder,  with  the  assumption  that  the 
plane  containing  the  source  and  the  axis  of  the  cylinder  is  perfectly  reflecting,  as  is  the 
surface  of  the  cylinder  itself.  The  technique  is  outlined  by  Huang  (1975),  extending  the 

plane-wave  results  derived  by  Friedlander  (1958).  More  difficult  problems  have  been 

• • 

attempted.  Uberall  et  al.  (1968),  McNicholas  et  al.  (1968),  and  Rudgers  and  Liberal  I (1970) 
derive  formulas  for  scattering  by  non-rigid  cylinders.  Carlson  and  Hung  (1974)  obtain 
some  results  onthe  difference  in  scattering  cross  sections  of  non-symmetric  profiles  as 
compared  with  that  of  the  cylinder. 

However,  we  shall  restrict  ourselves  to  the  more  idealized  problem  in  which  we 
con  derive  hard  quantitative  results.  The  analysis  proceeds  as  follows.  Let  all  distances 
be  made  non-dimensional  by  the  radius  of  the  cylinder  a,  time  by  q/c.  Let  Y be  the 
non-dimensional  distance  along  the  axis  of  the  cylinder,  and  D the  distance  of  the  source 
from  the  origin  on  the  axis.  The  distance  R between  the  source  and  the  microphone  is  thus 
given  by 

R^  = (X  + D)^  + Y^. 

This  geometry  is  diagrammed  in  Figure  6.2. 

Huang  (1975),  using  the  line-sourge  solution  of  Frienlander  (1958),  derives 
the  following  approximations  for  an  incident  unit  spherical  pressure  pulse.  The  line- 
source  (i.e. , two-dimensional  solution  is 


1 


F(X,f  - t ) = 0.924  [ - 1 KD^  - 1)  e (t  - Of* 

° “ (6.1) 

. H(t-t  )exp[-  0.280ej^'^^(t-O"^] 

Q C Q 

where  !■  is  the  fwo-dimensional  arrival  time  aX/c,  and  0 is  the  shadow  boundary, 
a c 

a — “11  “11  (A  oS 

0 ^ = It  “ COS  -g  “ COS  ^ . (0. 2) 

As  in  Section  5,  the  assumption  has  been  made  that  ^ '*  small,  thot  is,  that  the 

signal  has  not  long  passed  the  position  of  the  microphone  (X,Y).  The  solution  in  three 


dimensions  is  then  given  by 


A 1 F { [ t2“  y^“  t (1  + a)^(t^“  t^ ) ^ - t } 


1 J o il*  • I L -Q  g' 

P(X,Y,  t “ t ) =-  t(r  - r y J —5 5-^ 2 *** 

' ' ' a It  a **  , .2  .,2* 


“1  t^“i(l  +af  (t^-O) 


This  ejqjression  has  been  evaluated  by  numerical  computation,  and  the  various 
results  are  plotted  in  Figures  6.3,  6.4,  6.5,  and  6.6.  We  note  first  that  the  larger  the 
distance  from  the  cylinder  of  source  and  microphone,  the  more  the  pressure  trace  resembles 
a step  function,  that  is  the  smaller  the  influence  of  the  obstacle.  For  source  and  microphone 
close  in  to  their  respective  sides  of  the  hill,  the  diffraction  is  very  marked,  in  the  usual 
pattern:  the  amplitude  rises  slowly  from  zero  to  a much-delayed  maximum.  Which  point 
in  time  of  such  a signature  the  microphone  would  record  as  "arrival"  is  a matter  for  further 
investigation;  also  whether  or  not  two  different  microphones  would  give  the  same  answer. 
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Figure  6.3  Diffraction  of  a line-source  pulse  by  a cylindrical  ridge, 
computed  from  the  analysis  by  Friedlander  (1958),  The  ordinate  is  dimensionless 
pressure,  the  abscissa  is  dimensionless  time  after  arrival.  Each  curve  is  labeled 
by  a number  pair:  the  distance  of  the  source  from  the  center  of  the  ridge  and  the 
distance  of  the  microphone  from  the  center  of  the  ridge.  The  result  is  symmetric 
in  these  tv/o  quantities,  so  that  they  need  not  be  designated. 
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The  Simulation  of  a Howitzer  Blast 
Traversing  a 500-foot  Hill 

This  computation — the  first  of  its  kind  ever  attempted — is  structured  os  follows. 

The  spherically  symmetric  Howitzer  blast  is  computed  out  to  a radial  distance  somewhat 
less  than  500  feet.  At  this  point  the  data  from  this  computation  serves  as  input  into  the 
boundary  of  a two  dimensional  grid.  This  input  will  have  the  form  of  an  N-wove,  os 
represented  in  the  appropriate  profiles  of  Figure  7. 1.  The  two-dimensional  zones  are 
about  one  foot  on  a side,  so  that  smoothing  of  the  spherically  symmetric  solution  is 
necessary.  This  operation  wipes  out  the  spurious  numerical  oscillations  discussed  in 
Section  7,  so  that  they  do  not  appear  in  the  results  of  this  Section.  The  smoothing  is 

y 

done  occording  to  a scheme  that  conserves  both  energy  and  momentum. 

Since  the  solution  is  now  two-dimensional,  the  velocity  will  hove  two  components. 
A vector  velocity  plot  thereby  becomes  a meaningful  visualization  of  the  motion  patterns. 
These  must  be  viewed  with  the  foct  in  mind  that  the  vectors  represent  values  at  their 
base  points,  and  therefore  may  appear  on  occasion  to  penetrate  the  boundary.  The 
mathematical  condition  at  the  hill  is  zero  normal  velocity  component,  that  is,  perfect 
reflection. 

The  first  plot  of  Fig.  6.7  exhibits  the  grid  scheme.  In  this  plot  it  is  resting 
at  the  bottom  of  the  hill  on  the  side  of  the  source;  the  slope  of  the  hill  can  be  seen 
extending  up  to  110  feet  on  the  right.  As  the  pulse  progresses,  the  grid  is  moved  along 
at  a speed  that  keeps  it  ahead  of  the  front,  but  with  enough  zones  to  follow  what  form 
the  signal  will  take  at  any  given  point. 


The  plof  for  Cycle  50  shows  fhe  beginning  of  pulse  reflecHon  from  t^he  hill, 
particularly  pronounced  where  the  incidence  is  normal.  Cycle  50  shows  some  interesting 
features.  The  pattern  of  both  incident  and  reflected  waves  is  evident,  together  with 
their  complex  interoction.  But  also  in  this  cycle,  at  the  lower  left,  we  see  that  the 
reflected  wave  has  intersected,  and  been  re-reflected  from  the  grid  boundary.  Thus 
there  is  a triangle  in  this  portion  of  the  grid  that  contains  non-physical,  numerical 
results.  This  is  true  of  all  subsequent  cycles;  to  perform  the  computations  in  any  other 
way  vA>uld  be  prohibitively  expensive.  However,  this  region  of  physical  non -significance 
will  move  with  the  local  speed  of  sound  and  therefore  cannot  contaminate  more  than  a 
portion  of  the  grid.  In  Cycle  450,  for  example,  fhe  reinforcement  of  the  incident  wave 
by  tfie  reflected,  observable  olong  a line  from  the  wave  front,  is  physically  real.  At 
any  given  point  in  space  where  this  reinforcement  is  td<ing  place,  it  would  appear  os 
an  echo  if  there  had  been  a sufficient  period  of  silence  since  the  passage  of  the  front, 
or  os  a rumble,  if  not.  By  Cycle  850,  this  pattern  still  persists,  but  has  become  con- 
siderably more  complicated.  At  this  time  the  front  is  approaching  the  siwdow  limit, 
which,  as  calculated  above,  is  on  the  hillside  at  an  elevation  of  433  feet. 

By  Cycle  1100,  the  pulse  front  is  definitely  into  the  shadow,  and  the  process 
of  diffraction  has  begun,  although  scarcely  noticeable  at  this  stage.  Cycle  1300  shows 
a more  advanced  stage  af  diffraction:  the  amplitudes  are  overall  smaller,  and  they  rise 
from  zero  amplitude  at  the  front  relatively  slowly  to  their  maxima. 

This  pattern  is  still  more  in  evidence  at  Cycle  1909,  at  which  the  pulse  front  has 
just  reached  the  top  of  the  hill.  The  scale  of  fhe  velocity  vectors  of  the  earlier  plots 
has  been  retained,  and  the  diffracted  field  is  smooth,  regular,  and  slowly  varying. 
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In  Fig.  6.9  w«  see  a time  history  at  a point  just  five  feet  above  the  surface  on  top  the 
hill.  Here  the  pressure  trace  is  an  N-v/ave,  but  because  of  the  diffraction  it  is 
difficult  to  specify  a precise  arrival  time.  Even  at  this  distance  the  effect  of  non- 
linearity is  significant.  The  distance  from  the  source  is 

[1000^  + 500.5^  = 1118  ft 

The  arrival  time  accroding  to  linear  theory  for  the  sound  speed  1 1 16.36  ft/sec  would  be 
1. 002  sec.  We  see  from  Fig.  6. 9 that  it  is  closer  to  0, 96 

By  .Cycle  2850,  the  pulse  has  proceeded  about  100  feet  down  the  far  side  of  the 
hill.  The  diffraction  effect  by  Hiis  stage  is  very  significant.  The  pattern  is  such  that 
arrival  time  is  very  difficult  to  define.  The  scale  has  been  doubled  from  the  previous  plots, 
but  the  first  maximum  is  very  small  indeed.  The  signal  then  falls  off  to  negligiblity  and 
rises  to  a second,  i-einforced  maximum. 

In  the  plots  for  Cycles  3150  and  3650  the  scale  has  been  increased  by  a factor 
of  ten.  These  sliow  the  pulse  before  and  after  passing  the  45  degree  point  on  the  hill 
(354  feet  elevation).  To  show  the  effects  of  diffraction  we  compare  Hie  signal  signatures 
at  this  point  (Figure  6.8)  with  the  symmetric  point  on  the  other  side  of  the  hill,  that  is 
the  side  receiving  the  direct  radiation  (Figure  6. 10). 

The  signal  arrives  in  Figure  6.8  as  a true  pulse  with  a sudden  rise  and  the  beginning 
of  a decline  into  the  N-wave  configuration.  However,  by  this  time  of  passage  of  the  first 
pulse,  the  reflections  have  begun  to  arrive.  These  take  the  form  of  oscillations  with  a 
frequency  of  about  100  Hz.  The  signature  suggests  a lower  Frequency  component  olso,  but 
is  terminated  for  computational  reasons,  to  wit,  Hie  end  of  the  moving  grid.  It  should  be 


I 
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empfiasized  here  that  these  oscillations  are  physical  and  not  computational.  We  also 
note  that  pressure  and  velocity  are  precisely  in  phase,  a characteristic  of  the  direct 
linear  sound  propagation. 

Figure  6.10  presents  a very  great  contrast.  The  distances  of  the  two  points 
from  the  source  are  737  ft  and  1400  ft,  which  >vouId  result  in  a direct  radiation  amplitude 
reduction  of  about  47  per  cent.  The  actual  computed  reduction  of  the  amplitude  of  the 
first  pulse  maximum  is  seen  to  be  about  83  per  cent,  clearly  an  effect  of  the  shadow. 

The  structure  of  the  signal  has  been  totally  altered.  The  gradual  rise  from  zero 
to  the  first  maximum,  typical  of  the  diffracted  pulse,  mokes  impossible  the  field  detection 
of  true  arrival  time.  None  of  the  reflection  oscillations  are  evident;  and  the  second 
maximum  exceeds  the  first. 
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Figure  6.7  (following  pages).  Vector  plots  of  diffraction  of  Howitzer 
blast  by  a ridge. 


The  ridge  of  circular  cross-section  with  radius  of  500  feet,  is  centered 
1000  feet  from  the  howitzer.  Plotted  are  the  velocity  vectors,  scaled  according  to 
the  legend.  The  distance  coordinate  (in  feet)  is  X,  the  height  coordinate,  Y. 

Each  plot  shows  the  pulse  front  and  some  distance  behind  it.  The  point  X = 750, 

Y = 433  on  the  hill  begins  the  shodow  zone.  See  text  for  further  details. 
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7.  THE  PROPAGATION  OF  A HOWITZER  BLAST 

The  serious  theoretical  study  of  pulse  propagation  from  an  intense  explosion 
begain  as  classified  research  during  the  second  World  War.  The  pioneering  discovery 
of  self-similar  solutions  to  the  spherical  problem  was  the  independent  work  of  Taylor 
(1946,  1950)  and  of  Sedov  (1959,  Ch.  4).  Building  upon  this  work,  Whithorn  (1950; 
1974,  Chs.  6,9)  developed  a technique — later  termed  "non linearization” — that 
calculates  approximations  based  on  large  distance  from  the  source,  but  tckes  into 
account  the  fact  that  the  true  characteristics  are  not  of  the  linearised  form 


f “ ~ constant,  but  are  some  function  T(R,t)  = const,  to  be  determined.  To  use 

the  linearized  form  might  give  good  results  near  the  source,  but  would  lead  to  serious 
error  for  large  R owing  to  the  divergence  of  the  true  characteristics.  However,  for 
explosions  of  the  magnitude  that  concern  us  here,  the  overpressure  orjd  the  particle 
velocities  are  small  when  compared  respectively  to  the  mean  pressure  and  the  sound 
speed,  except  very  near  the  source.  The  disturbance  is  thus  always  paradoxically 
propagating  like  a linear  distance  in  one  sense  and  like  a non-linear  one  in  another 


sense. 


The  theoretical  treatment  (Whitham,  1950,  pp,  579  ff. ) takes  as  a model 'an 
expanding  sphere,  initially  of  radius  d,  of  uniform  high  pressure  p®/p  » 1 where  p 
is  atmospheric.  The  uniformly  valid  approximation  for  large  distance  R is  then  os  follows. 
The  shock  front  has  the  representation 


c t = R - b(C7iR)^ 


2 p 2 

where  b =0,84  ^ d . The  constant  b cannot  really  be  measured,  since  the  expanding- 
sphere  model  is  not  an  adequate  representation  of  the  actual  blast.  However,  if  we 


assume  a pressure  of  say,  40  ofmospheres  in  a sphere  of  one  mefer  radius,  we  obfain 
b = 5. 8 m and  for  a distance  R of,  say  5 kilometers,  we  have 


b(enR)^  =16m  (7.2) 

For  a scihd  speed  c = 33O  m/sec,  this  would  represent  a correction  of  0.03  seconds  on 
o 

a linear  arrival  time  of  15. 15  seconds. 

The  reason  for  this  is  clear:  during  its  highly  nonlinear  expansion,  the  shock 
travels  with  a speed  in  excess  of  the  sound  speed.  The  equation  (7. 1)  shows  that  the 
linear  propagation  is  not  quite  effectively  achieved,  even  after  a considerable  distance 
and  time  of  propagation. 

However,  the  picture  is  still  more  complicated.  After  the  first  overpressure  was 

has  passed  from  the  explosive  gases,  a rarefaction  wave  is  propagated  back  into  them. 

This  is  followed  by  a new  second  shock,  a second  rorefaction  wove,  a third  shock,  ond 

so  on  until  the  explosive  gases  have  been  effectively  reduced  in  energy. 

Theory  tells  us  that  the  second  shock  propagates  as 

c t = R + b(0nR)^,  (7.3) 

o 

a counterpart  to  (7. 1).  However,  the  details  of  the  shock  formation  depend  essentially 
on  the  time-geometry  of  the  expanding  sphere,  and  this  is  precisely  the  most  unrealistic 
aspect  of  the  model.  Thus  one  proper  test  of  the  theoretical  results  is  a full-scale 
computation  of  the  explosion  and  of  the  subsequent  shock  propagation. 

The  Numerical  Simulation  of  a Spherical  Blast 
There  hove  been  a number  of  computer  simulations  of  spherical  blast.  One 


highly  successful  example  is  that  of  H.  L.  Brode  (1959).  This  particular  piece  of 


research  by  Brode  represents  the  culmination  of  a series  of  similar  efforts  by  himself 
and  others.  We  have  used  the  Afton  Code  (Trulio  et  al. , 1976)  to  simulate  a howitzer 
blast,  using  as  input  the  parameters  associated  with  the  WSMR  PASS  operation,  together 
with  the  most  recent  equations  of  state  for  explosives.  The  basic  assumption — discussed 


/ 


in  later  sections — is  that  dispiite  the  essentially  directiortal  nature  of  a howitzer  explosion, 
for  results  concerned  with  even  moderate  distances  from  the  source,  the  event  may  be 
idealized  as  spherically  symmetric.  This  simplifying  geometry  permits  us  to  compute  the 
expansions  and  rarefactions  with  a grid-network  sufficient  to  resolve  the  details.  In  the 
early  stages  of  the  explosion,  for  example,  it  was  necessary  to  use  a zone-dimension  of 
7.5  cml  The  physical  quantities  of  the  computation  were  provided  by  personnel  of  the 
Atmospheric  Sciences  Laboratory,  WSMR  (Blomer.r>,  personal  communication).  Together 
with  other  derived  parameters,  they  ore  as  follows. 


Initial  Conditions 


In  undisturbed  air: 
density 
pressure 
sound  speed 
In  the  high  explosive: 

radius  of  HE  sphere 
density 

specific  energy 


p^  = 1 . 225  kg/m^ 

p = 1.013  bars 

c = 340.25  m /s 
o 

Rg  = 0. 1634  m 

Pg  = 1557  kg/m^ 

eg  = 4. 2945  x 10^  joule  Ag 


(7.4) 


The  cocnputort^ion  obviously  requires  one  set  of  equations  for  the  explosion  and  the 
subsequent  movement  of  the  exploded  gases,  and  another  for  the  resulting  motions 
of  the  air.  Between  these  two  gaseous  systems  there  will  remain  an  interface. 

The  dynamic  equations  for  the  high  explosive  are  those  of  Taylor  (1950), 
with  the  exception  of  the  equation  of  state.  Because  the  constants,  and  even  the  form, 
of  the  equation  of  state  for  explosives  are  best  determined  empirically,  we  have  used 
for  this  purpose  the  recent  results  of  Lee  et  al.  (1968)  as  the  most  satisfactory  version 
to  date. 

The  air  was  treated  as  a gamma-low  gas, 

Ap  = y £ Ap  , 

and  integrated  according  f©  the  Afton  Code.  The  finite  difference  equations  of  this 
code  are  complex  but  fully  documented  by  Trulio  et  al.  (1976)  in  the  User’s  Manual. 
Two  runs  were  performed  for  the  spherically  symmetric  explosion: 

(1)  Using  special  zones  of  7.5  cm  length,  the  integration 
was  carried  out  to  1.5  seconds. 

(2)  Using  spacial  zones  of  15  cm,  the  integration  was  carried 
out  to  6.735  sec.  This  result  was  used  as  the  input  to  the 
computation  of  diffraction  by  a hill,  reported  in  Section  6. 

It  may  be  mentioned  at  this  point  that  to  our  knowledge^  these  computations  of  a blast 
wave  ore  here  carried  out  an  order  of  magnitude  further  in  time  (and  space)  than  any 
other  yet  performed.  The  computations  of  Brode  (1959),  for  example,  are  terminated 
(corresponding  to  the  data  for  our  problem)  at  about  0.12  sec. 
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is  still  active  dynamically  in  a very  complicated  pattern.  TTie  shock-producing  process 
may  or  may  not  still  be  in  operation,  but  the  magnitudes  are  no  longer  sufficient  to  be 
detected  above  the  computational  noise. 

This  noise — that  is,  truncation  error — is  evident  in  the  periodic  wiggles  back 
of  each  shock;  and  a word  concerning  this  problem  is  in  order.  The  shock  front,  having 
the  theoretical  form  of  a step  function,  cannot  be  resolved  In  its  propagation  by  finite 
differences.  The  cost  of  discretizing  the  computation  is  that  spurious  oscillations  are 
set  up,  and  the  amplitude  of  the  shock  energy  maximum  is  diminished  by  the  energy  of 
these  little  wavelets.  This  situation  has  been  examined  in  some  detail  by  Wurtele  (1961), 
who  presents  analytical  solutions  for  differential-difference  equations  that  enc4>le  us  to  see 
the  structure  of  the  numerically  introduced  error.  By  the  methods  of  this  paper,  we  con 
specify  that  at  a time  of  51  ms,  the  numerical  oscillation  should  have  a wave  length  of 
6. 1 zones.  The  dump  at  51.492  shows  that  this  is  indeed  the  wavelength.  However,  in 
order  to  check  further,  we  repeated  the  computation  with  a reduced  zone  length  of  7.5  cm. 
The  results  are  contained  in  Figure  7. 2.  A dump  at  the  same  time  (51  ms)  shows  a finer 
resolution  of  the  shocks,  to  be  sure,  but  witli  double  the  number  of  oscillations  per  unit 
distance.  This  demonstrates  beyond  question  that  the  oscillation  is  numerical  and  not 
physical. 

It  would  be  simple  to  remove  these  spurious  oscillations  by  smoothing — a running 
average  over  a dozen  zones,  for  example — but- this  would  be  a purely  cosmetic  operation, 
with  no  information  gained  and  with  new  problems  introduced  which  we  would  understand 
less  well  than  we  understand  the  existing  numerical  error.  For  this  reason,  we  hove 
graphed  the  computed  solution  with  no  modification  at  any  time. 
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The  se^  of  graphical  output  in  Figure  7. 1 show  the  simulation  of  the  blast. 

The  upper  figure  (labelled  stress)  on  each  page  is  the  overpressure,  the  middle  is  the  parcel 
displacement,  and  the  bottom  is  the  (radial)  velocity.  Each  graph  represents  the  distri- 
bution of  the  quantity  in  space  at  a given  time.  Care  must  be  taken  in  reading  the 
magnitudes,  since  both  the  ordinate  and  the  abscissa  scales  change  with  time  in  order 
that  a satisfactory  graph  be  presented.  The  zero-point  also  is  variably  placed. 

With  these  facts  in  mind,  we  may  examine  the  results.  The  first  shock  is 
evident  in  DUMP  1,  at  0.024  ms.  The  dumps  at  0.693  and  1. 121  ms  show  the  sub- 
sequent extreme  rarefaction.  In  the  latter  graph  the  pressure  appears  to  be  zero— 
i.e.  a negative  one-bar  overpressure — and  reference  to  the  digital  output  shows  it  to 
be  indeed  almost  but  not  quite  a vacuum. 

Recallir^  that  this  is  a spherical  rather  than  a linear  blast-wave,  we  note  that 
the  rarefaction  cannot  be  propagated  backward  beyond  the  origin,  and,  according  to 
the  above-mentioned  theory,  a second  shock  must  follow.  This  is  first  evident  at  - 
6.313  ms:  the  first  shock  has  propagated  to  about  22  feet  from  the  source;  and  the 
rarefaction  extends  from  about  15  feet  back  to  5 feet,  where  the  second  shock  begins. 

This  process  is  repeated  twice  more.  Reference  to  the  dump  at  51 .492  ms  shows  an 
interesting  situation.  Four  shocks  are  evident,  at  82,  55,  46,  and  36  feet.  Each 
shock  is  smaller  than  its  predecessor.  The  first  has  still  a significant,  though  greatly 
reduced,  overpressure  (0.07  bars);  and  the  subsequent  shocks  exist  within  a general- 
rarefaction  region,  with  each  sub-peak  pressure  just  zero  or  slightly  positive.  At  about 
seven  feet  one  notes  the  interface  of  the  high  explosive  and  the  gommo-law  air.  The  HE 
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We  may  now  look  at  a later  time,  say  176.507  ms.  At  this  point  the 
computation  within  the  HE  was  abandoned  os  no  longer  contributing  significantly  to 
the  region  of  interest  behind  the  shock.  The  first  shock  itself  has  reached  225  feet  and 
its  maximum  overpressure  has  fallen  to  0.05  bars,  or  50  mb.  The  second  and  third 
shocks  are  still  detectcble  above  the  computational  noise,  but  the  fourth  is  now  lost. 
The  velocity  maximum  is  about  13  m sec  ^ ^ small  compared  to  the  speed  of  sound, 
but  certainly  large  for  a sound  wave. 

The  lost  dump,  at  6701 . 280  sec,  shows  the  shock  front  at  7520  ft.  Only  two 
shocks  are  now  distinguishable.  The  spurious  oscillation  following  the  first  shock  is 
about  40  zones  in  wave  length.  Overpressure  is  about  one-half  millibar  and  velocity 
15  cm/sec.  It  is  this  state  of  the  system  that  served  os  input  to  the  computation  of 
scattering  by  a hill  in  Section  6. 

In  order  to  make  connection  with  theory  let  us  consider  a relatively  late 
dump,  soy,  70.  This  is  at  time  2.766  seconds,  using  the  sound  speed  (7.4),  we  hove 

c t = 941. 13  m. 
o 

Now  owing  to  numerical  dispersion,  we  have  an  interval  bounding  the  arrival  time, 
rather  than  a precise  time  of  arrival.  An  understanding  of  the  numerical  effects  would 
lead  us  to  place  the  shock  front  at  DUMP  70  at  3130  feet  or  954. 02  m.  It  could  not 
be  placed  at  less  than  3115  feet  or  949.45  m.  These  distances  translate  to  time  dif- 
ferences of  0.03788  and  0.02445  sec  respectively.  Thus  in  the  simulation,  os  in  the. 
nort-lineor  theory,  the  signol  arrives  meosurably  ohead  of  the  orrivol  time  predicted 
by  the  bosic  sound-ranging  assumptions. 
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If  we  wish  to  establish  even  more  detailed  correspondence  between  simulation 
and  analysis,  we  may  apply  (7. 1)  to  the  simulation  results.  Using  R = 954.02  m and 
t = 2.766  sec,  we  obtain  a value  for  b, 

R - c t 

b=  =6. 4m 

(0nRf 

This  is  not  far  from  our  tentative  value  of  b = 518  m,  calculated  on  the  basis  of  a pressure 
of  40  atmospheres  in  a one-meter-radius  sphere.  If  one  examines  the  early  dumps,  he 
will  see  that  although  the  expanding -sphere  model  is  not  really  applicable,  if  one  were 
forced  to  use  it,  the  values  assumed  would  be  defensible. 

At  any  rate,  we  see  that  the  arrival  time  is  distorted  by  the  essential  non- 
linear! ly  of  the  shock  propagation.  It  remains  to  determine  how  great  an  error  is 
introduced  into  sound-ranging  results  by  this  distortion. 

Implications  for  Sound  Ranging 

Using  some  rather  gross  assumptions,  we  may  estimate  the  ranging  error  involved 
with  the  nonlinearity.  Let  us  use  (7. 1)  v/ith  b = 6.  Then  if  we  apply  (7. 1)  to  the 
arrivals  at  R^  and  R^,  and  subtract  the  results,  we  hove 

c^At  = AR  - b 3^-  [^(R^)]^} 


(7.6) 


where  At  = t^  - t^ , and  AR  = Rj  ■ R^j.  We  now  use  the  approximations — solely  for  the 
purpose  of  the  present  discussion  — 

AR 

072  (1  + -K—  ) « 072  R, 

R]  1 


and 


AR 

p-  « 1 . 


r 


Then  (7. 6)  becomes 


2n(l+^) 

V 


= AR  n X-  ■ ] 

2R,(enR^r 

TTiis  applies  equally  to  the  signal  arrivals  at  R and  R . Hence  we  may  write 

I o 


s.  = (1  + e ) ff. 


1=2,3 


where,  as  before,  a. -c  At.,  and 
I o I 

1+€=[1 T^"' 

2R^(C7iRp^ 

_3 

It  is  easily  estimated  that  for  b - 6 m,  the  value  of  e will  be  about  10  . If  we 

substitute  (7.7)  into  the  ranging  formula  (2.3a),  we  obtain 


2R  = 


2d2-((r32-2a2^)(l+e)^ 
(1+  c)(ff3  - 2a ^ 


or  to  the  first  order  in  € 


R - R (1  -c)  - € 


3-^7 

■’3-^2 


The  correction  of  one-tenth  per  cent  on  R^  is  perhaps  bordering  on  negligibility, 
but  the  additive  term  would  normally  be  a correction  three  or  more  times  as  large. 

It  would  seem  that  a correction  factor  as  in  (7.9)  should  be  built  into  the  ranging 
formulas  and  computational  code.  The  approximations  made  in  deriving  (7.8)  might 
affect  the  result,  but  (7.9)  would  always  lead  to  a more  accurate  ranging  estimate. 
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We  are  not  quite  in  a position  to  recommend  a final  formulation  for  this 
correction.  The  evidence  seems  to  be  that  a value  of  six  for  b is  adequate,  but  further 
computation  would  be  desiroble. 

The  value  of  in  (7.8)  is,  of  course,  unknown.  In  theory,  an  estimate  could 
be  made,  computed,  and  then  an  iteration  procedure  performed.  However,  in  view 
of  the  smallness  of  the  correction  in  practice,  this  seems  not  worthwhile.  The  first 
estimate  is  probably  adequate. 

I 

The  amplitude,  as  well  as  the  arrival  time,  is  always  measurably  nonlinear 
in  behavior.  For  example,  the  peak  overpressure  at  any  distance  R is  given  by 


P 

P 


Ad 


(7. 10) 


o R(  Cn  R)* 

where  A and  B are  to  be  empirically  determined  Miles  (1967)  has  gathered  the  available 
data,  and  for  comparison  with  linear  theory  determined  an  "apparent  exponent  of  decoy" 
q such  that 

P~R“^.  (7.11) 

In  lineor  theory,  of  course,  q - 1.  Alrhougli  behavicr  according  to  (7. 10)  can  never 
satisfy  a power  law,  it  turns  out  that  for  typical  sound-ranging  distances  R^  the 
value  of  q is  relatively  constant,  about  1.07.  This  and  the  other  empirically-determined 
values  could  be  introduced  into  (7. 11)  so  that  we  can  replace  (7.8)  by 


c = 


K 

.1.07 


(7.12) 


expressed  in  meters,  and  the  constant  K is  to  be  determined  by  field  experiments. 
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Figure  7. 1 The  Numerical  Simulation  of  a Howitzer  Blast. 

Each  page  shows  special  profiles  of  three  quantities  at  a given  time. 

The  top  profile  is  that  of  the  stress,  that  is,  the  total  pressure.  The  second  is  the 
parcel  displacement,  and  the  bottom  is  the  velocity.  Since  the  explosion  is  assumed 
sphericals  the  velocity  is  positive  outward,  negative  inward,  with  the  origin  at  the 
center  of  the  explosion.  The  distance  from  this  center  is  the  abscissa,  labeled  "radius". 

The  scales  and  the  zero  positions  are  chosen  to  obtain  a readable  graph 
and  should  be  consulted  in  viewing  each  profile. 

The  time  of  eoch  page  is  given  at  the  top  right  in  milliseconds. 
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Figure  7.2 

Same  as  Figure  7. 1,  but  with  space  zones  of  7, 5 cm,  that  is,  twice  the 
resolution  of  the  previous  figure. 
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8.  ANALYSIS  OF  PASS  DATA 


Twenty  separate  Howitzer  firing  events  were  studied  from  the  data,  for  two 
of  which  the  actual  microphone  traces  were  available.  These  events  were  chosen 
by  referring  to  the  wind  and  temperature  profiles  from  the  500  foot  instrumented 
tower  located  near  the  two  Howitzers.  Conditions  of  near  calm  were  preferred  for 
all  events,  except  for  the  eight  on  12-5-74,  so  that  wind  velocities  would  not 
significantly  enter  into  the  range  calculations.  The  eight  events  on  12-5-74  were 
chosen  as  the  highest  surface  wind  speeds  available  to  test  the  wind  correction 
calculation.  Examples  of  the  three  typical  vertical  profiles  of  temperature  (lapse, 
isothermal  and  inversion)  were  included  in  the  twenty  events  studied. 

We  have  tabulated  each  event  studied  in  a format  known  as  the  Event 
Summary.  The  uppermost  section  of  each  event  summary  identifies  the  three  micro- 
phones chosen  for  the  baseline,  the  day,  the  event  number,  the  gun  number,  the 
time  of  Howitzer  firing,  expressed  as  Mountain  Standard  Time,  to  the  nearest 
millisecond  and  the  distances  between  microphones.  The  distance  between  the 
closest  and  the  middle  microphone  is  labelled  d^,  that  between  the  farthest  and 
the  middle  microphone  d , and  that  between  the  closest  and  farthest  microphone 

is  D — d + d . When  d and  d are  approximately  equal,  their  numerical  average 
2 2 o 

is  expressed  simply  as  d and  the  total  distance  D is  not  needed.  The  next  section 
down  on  each  event  summary  includes  the  speed  and  direction  of  the  wind  W at 
25  feet  elevation,  the  air  temperature  T at  26  feet  elevation,  the  calculated  speed 
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of  sound  C,  ttie  moch  number  of  fhe  wind  M,  and  the  angle  9 between  the  microphone 
baseline  and  the  wind  vector.  This  angle  is  measured  from  the  for  end  of  the  baseline 
back  in  tov«ird  the  Howihter  to  the  head  of  the  wind  vector.  The  third  section  in  each 
event  summary  is  a self-explanatory  table  directly  comparing  the  raw  microphone  time 
data  to  the  predicted  or  correct  time  information  using  the  temperature  calculated  speed 
of  sound.  The  measured  distance  between  Howitzer  number  one  and  any  microphone  is 
given  as  rj  and  the  measured  distance  between  Howitzer  number  two  and  any  microphone 
is  given  as  rj.  The  raw-data  microphone  times  are  given  as  t.  The  time  error  column  shows 
a negative  number  for  late  arrivals  (recorded  later  than  calculated)  and  a positive  number 
for  early  arrivals.  The  bottom  section  of  each  event  summary  compares  true  time  differences 
between  microphones  with  raw  data  time  differences,  compares  measured  values  of  and 
S (see  wind  correction  and  ranging  diagrams)  with  values  for  S and  S calculated  from 

it  W 

raw  time  data  and  the  temperature  calculated  speed  of  sound,  compares  the  row-data- 

calculated  angle  between  ;he  microphone  base  line  and  line  between  the  farthest  microphone 

and  the  Howitzer,  Oi-,  with  the  wind-corrected  calculation  of  that  angle,  di  , , and  the 
^3  3wind 

true  angle,  True  ij»-,  and  compares  the  row-data-calculated  distance  from  nearest  microphone 


to  Howitzer,  R , with  the  wind-corrected  calculation  of  that  distance,  R . ,. 

o wind 

On  the  adjoining  two  pages  to  this  written  summary  we  have  compiled  the 
results  of  the  twenty  events  studied.  One  quick  glance  reveals  the  rather  serious 
problem  of  uncertain  and  missing  data.  This  category  accounts  for  nearly  half 
of  all  possible  microphone  reports  in  the  study.  The  criterion  for  deciding  which 
arrivals  were  early,  late  and  correct  was  very  simple.  I considered  as  correct 
results,  those  that  exactly  (to  the  nearest  millisecond)  corresponded  to  the  predicted 
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or  true  values  obtained  by  using  the  temperature  calculated  speed  of  sound. 

Under  these  circumstances  one  would  not  expect  there  to  be  many  "correct" 
reports;  indeed  only  one  was  found.  But  a good  microphone  might  be  one  that 
has  an  even  distribution  between  late  and  early  arrivals,  microphone  B3  for 
example.  After  one  discards  the  missing  data,  only  278  good  resclts  out  of  a 
possible  520  are  seen.  Of  these,  a striking  72.3^  of  the  reports  are  too  CarlVt 
This  compilation  of  twenty  events  would  seem  to  filter  out  any  systematic  time 
error  due  to  wind  direction  because  nearly  all  wind  directions  are  included  in 
these  results.  Another  possibility  for  systematic  error  would  be  an  error  in  the 
' timing  mechanism  of  the  Howitzer  firings  leading  to  a preponderance  of  "early" 

reports.  This  type  of  error  should  not  affect  the  sound  ranging  calculations  since 
the  absolute  error  in  the  arrival  times  is  not  important  in  itself.  It  is  the  time 
differences  between  microphones  that  is  important  and  this  is  not  affected  in  any 
way  by  an  error  in  the  determination  of  the  absolute  time  of  Howitzer  firing  t^. 

Results  of  the  Survey 

When  one  surveys  the  event  summaries,  he  cannot  fail  to  be  impressed 
with  the  very  large  errors  sometimes  produced  by  the  ranging  formulas.  In  not  a 
few  events,  the  error  in  the  range  was  of  the  same  size  as  the  range  itself.  These 
cases  stand  in  stork  contrast  to  those  in  which  the  formulas  give  highly  satisfactory  j 

results,  quite  sufficiently  accurate  for  the  purpose  for  which  sound  ranging  is 
intended. 
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We  hove  made  every  effort'  to  identify  the  cause  of  error  in  those  cases 
^ in  which  the  error  was  significant.  A brief  summary  of  these  efforts  follows. 

(1)  Is  it  possible  that  a major  source  of  error  could  be  the  ambiguity  of 
the  microphone  traces  in  indicating  arrival  times,  or  in  the  misreading  of  these? 
Here  we  were  able  to  investigate  two  of  the  20  events  for  which  we  had  been 
provided  with  copies  of  the  microphone  traces  for  the  microphones  D4,  D5,  A3,  A6, 
A7,  B2,  and  D7.  The  two  events  were  on  Day  323  (101^4,  events  13  and  34j, 

It  will  be  noted  that  range  errors  were  as  large  as  88  per  cent  in  these  events. 
However,  it  is  not  possible  to  ascribe  any  appreciable  error  to  ambiguities  or 
misreadings  of  arrival  times  from  the  traces.  In  every  trace  the  arrival  time  was 
clearly  indicated,  and  had  been  properly  read.  From  this  tiny  sample  we  cannot 
conclude  that  the  traces  were  never  a source  of  error;  but  we  can  see  that  a well- 
marked  and  well-read  trace  can  give  highly  erroneous  results. 

(2)  Wind  correction.  As  stated  above,  we  had  selected  out  the  cases  of 
no  wind  or  very  light  wind,  so  thot  at  the  outset  we  can  state  that  large  errors  can 
result  inder  almost  completely  calm  conditions.  If  we  now  turn  to  consideration  of 
our  cases  selected  on  the  basis  of  relatively  high  wind,  we  can  survey  the  results  of 
applying  the  wind  correction.  These  are  not  encouraging.  In  roughly  half  of  the 
cases,  the  wind  correction  increased  the  ranging  error;  in  no  case  did  it  turn  a 

'poor  result  into  a good  result.  Tbe  Mach  numbers  in  all  of  these  examples  were 
quite  small  (less  than  0.04);  fhe  PASS  operation  was  conducted  primarily  under 
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conditions  of  little  weather.  When  the  wind  was  high,  a large  percentage  of 
the  microphones  were  missing  or  indicated  as  uncertain.  In  the  remaining  cases, 
those  examined,  we  conclude  that  the  wind  represents  a small  signal  in  the  general 
noise  level  of  whatever  processes  are  entering  into  the  error.  In  such  conditions, 
a wind  correction  would  not  be  expected  to  produce  improved  results. 

(3)  Terrain.  Reference  to  the  topographic  maps  of  the  firing  area  indicates 
no  physical  features  suggesting  that  diffraction  by  terrain  could  play  any  role  in  the 
errors  recorded  in  the  Events  Summary. 

(4)  Non-linear  and  directional  blast  effects.  Here  we  refer  to  Section  7 
in  particular,  formulas  (7.7)  and  (7. 12): 

e = and  S.  = (1  + e)  a. 

The  theory,  together  with  some  approximations,  suggests  that  the  constant  K has 
a value  in  the  neighborhood  of  1 to  5.  However,  this  assumes  a perfect  spherical 
wave,  which  may  not  be  an  adequate  model  for  the  non-linear  effect  of  a Howitzer 
blast.  Accordingly  we  have  followed  a common  procedure  in  selecting  a formula 
such  as  the  above  that  has  theoretical  grounding,  but  allowing  experments  to 
determine  the  value  of  the  constant.  For  this  experiment  we  selected  the  A 
microphones  as  being  most  directly  in  the  line  of  Howitzer  fire.  Reference  to  the 
Event  Summary  will  show  that  for  these  cases,  a value  of  K between  10  and  50, 
but  closer  to  10,  gave  improvement  in  the  ranging  estimates.  In  some  cases,  this 


r 


improvemenf  was  marked.  This  does  not  guarantee  that  similar  improvement 
can  be  made  on  independent  data  using  this  range  of  values  for  K,  but  it  does 
suggest  that  such  a correction  is  well  worth  trying  if  such  data  become  available. 


When  all  of  the  above  factors  have  been  taken  into  consideration, 
a large  residue  of  cases  of  large  error  remains.  When  (i)  the  signal  arrival  is 
inexplicably  early  at  the  first  microphone  and  (ii)  the  signal  sufficiently  irregular 
in  its  arrival  at  the  other  microphones  to  produce  an  error  of  100  per  cent,  then 
we  are  compelled  to  ascribe  the  errors  to  the  timing  instrumentation.  No  phenomenon 
of  Hie  atmosphere  or  feature  of  the  terrain  could  account  for  this  pattern  of  results, 
which  was  observable  in  a surprisingly  large  prop-.'rtion  of  the  events  examined. 
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Day 


Events  Studied 


Event  Number 


11-8-74 

(312) 

12 

11-8-74 

(312) 

14 

11-8-74 

(312) 

92 

11-15-74(319) 

95 

11-18-74(322) 

41 

11-19-74(323) 

13 

11-19-74(323) 

34 

12-2-74 

(336) 

9 

12-2-74 

(336) 

94 

12-2-74 

(336) 

176 

12-3-74 

(337) 

91 

12-5-74 

(339) 

21 

12-5-74 

(339) 

33 

12-5-74 

(339) 

69 

12-5-74 

(339) 

83 

12-5-74 

(339) 

85 

12-5-74 

(339) 

87 

12-5-74 

(339) 

97 

12-5-74 

(339) 

99 

12-7-74 

(341) 

230 

Twenty  Separate  Events  Studied 


Microphone 

Arrival  Times 
Too  Short 

Arrival  Times 
Too  Long 

Arrivals 

Correct 

Uncertain  or 
Missing  Data 

(5)  N5 

9 

3 

0 

8 

(6)  N6 

4 

0 

0 

16 

(7)  N7 

3 

0 

0 

17 

(8)  N8,  El 

12 

1 

0 

7 

(11)  E2 

11 

3 

0 

6 

(12)  Al,  W3 

6 

5 

0 

9 

(13)  A2 

10 

3 

0 

7 

(14)  A3,  D3 

10 

2 

0 

8 

(15)  A4 

4 

0 

0 

16 

(16)  A5 

1 

0 

0 

19 

(17)  A6 

12 

0 

0 

8 

(18)  A7 

14 

2 

0 

4 

(19)  A8,  E3 

13 

1 

0 

6 

(21)  D4 

12 

0 

0 

8 

(24)  D5 

3 

8 

0 

9 

(26)  Bl,  W6 

7 

3 

0 

10 

(27)  B2 

8 

4 

1 

7 

(28)  B3 

6 

6 

0 

8 

(30)  B5 

3 

5 

0 

12 

(31)  B6,  D6 

8 

5 

0 

7 

(32)  B7 

13 

4 

0 

3 

(33)  B8,  E6 

11 

6 

0 

3 

(34)  W7 

6 

2 

0 

12 

(36)  E7 

3 

5 

0 

12 

(37)  SI,  W8 

7 

3 

0 

10 

(44)  S8,  D8, 

E8  5 

5 

0 

10 

Totals 

^ based  on  all 
520  possibilities 

^ based  on  only 
good  reports 
(278) 


201 

76 

1 

(38. 7^) 

(14.6^) 

(0.^) 

(72. 3^) 

(27.3^) 

(0.4f6) 

( 
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Day:  11/8/74(312)  B6  (31) 

Event:  12  ^32) 

Gun:  2 ^^3) 


t = 
o 

14:00:0.335  MST 

d = 4,430.30* 

W = 

2.6  kts  (4.39  ft/sec) 
from  330°.  2 at  25'  elev. 

T = 52.5°F  (11.4°C)  at  26' 
C = 1,109.9  ft/sec 

elev. 

M = 

0. 0039565 

0 = 

323°. 0 

Mike 

True  Distance 
(feet) 

True  Time 
r^c  (sec) 

Time  Data 
t(sec) 

Time  Error 
r/c-t  (sec) 

B6 

9,405.76 

8.474 

8.476 

-.002 

B7 

12,915.59 

11»637 

11.638 

-.001 

B8 

16,864.62 

15. 195 

15.208 

-.013 

True  (^2"*^]^  ~ ^^2  sec 

^^2"^1^  Data  = 3. 162  sec  (no  error) 


True  ~ 


6.720  sec 


(t^-tj)  Data  = 6.732  sec  (long) 
True  $2  = 3,509.83' 

TrueSg  = 7,458.86' 

DotaS^  = 3,509.50'  \ 

DataSg  = 7,471.85*  ) 

Wind  $2  =3,523 
Wind  = 7,499 


= 8,898.92' 


.50'  ) 
.85'  i 


wind 


= 8,653.50' 


True  4^  - 

3;irs' 


-*■  5.39^  low 
-►  8. 0O»  low 
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Day: 

11/8/74  (312) 

W6  (26) 

Event 

‘:  12 

W7  (34) 

Gun: 

2 

W8  (37) 

t = 
o 

14:00:  0.335  MST 

d = 4,430.33' 

W = 

2.6  kts  (4.39  ft/sec) 

T = 52.5°F  (11. 

4°C)  at  26' 

elev. 

from  330°.  2 at  25'  elv. 

C = 1,109.9  ft/ 

sec 

M = 

0. 0039565 

0 = 

o 

• 

o 

True  Distance 

True  Time  Ti 

ime  Data 

Time  Error 

Mike 

(feet) 

r^c  (sec) 

t(soc) 

r^/c-t  (sec) 

W6 

16,898.89 

15.226 

15.189 

.037 

W7 

19,087.77 

17.198 

17.164 

.034 

W8 

21,962.93 

19.788 

19.760 

.028 

True  - 1.972  sec 

Dota  = 1.975  sec  (long) 

True  “ 4.563  sec 

Data  = 4.571  sec  (long) 

TrveSj  = 2,188.88'  - 

TrueSg  = 5,065.04' 

Data  $2  = 2,192.05'  ) ^ 

Data  S = 5,073.35'  ( ® 

nS 


Wind  $2  ~ 
Wind  S.  = 


= 2,202.60'  ) 
= 5,094.45'  ) 


R . 16,688.52' 

wind 


Twt  V|)^=  ifs-.o 

¥K7 


= 16,776.88'  ^ 0.722^  low 


.•>  1 . 24-^  low 
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Day:  1 1/8^4  (312) 

Event:  14 

Gun:  2 

t = 14:02:05.367  MST 


W = 2.6  kts  ^.39  ft/sec) 

from  330  . 2 at  25'  elev. 

M = 0.0039565 
0 = 323°.  0 


B6  (31) 

B7  (32) 

B8  (33) 
d = 4,430.30' 


T = 52.5°F  (11.4°C)  at  26'  elev. 


C = 1,109.9  ft/sec 


Mike 

True  Distance 

Tj  (fe»0 

True  Time 
(sec) 

Time  Data 
t(sec) 

Time  Err 
( 

1.) 

B6 

9,405.76 

8.474 

8.545 

-.071 

2.) 

B7 

12,915.59 

11.637 

11.633 

+.004 

3.) 

B8 

16,864.62 

15.195 

15.202 

-.007 

True  (t^-t^)  = 3. 162  sec 

^^2~^1^  Data  = 3.088  sec  (short) 

True(t2-t^)  = 6.720  sec 

(ta-f^)  Data  = 6-657  sec  (short) 

True  $2  = 3,^09.83* 

True  S_  = 7,458.86 

w 


Tv-ot 

Y,=  :i5:s- 


Data  $2  ~ 3,427.37' 
DataS  = 7,388.60' 
Wind  $2  = 3,441.37') 
WindSa=  7,416.60'! 


= 7,640.05' 


R = 7,431.93’ 

wind 


18.8^  low 
2l.Dib  low 


Day:  11/8/74(312)  W6  (26) 

Event:  14  W7  (34) 

Gun:  2 W8  (37) 


• = 14:02:05.367  MST  d = 4,430.33' 

o 


W = 

2.6  kts  (4.39  ft/sec) 

T = 52 

'.5°F  (11.4°C)  at  26' 

elev. 

from  330®.  2 at  25'  elev. 

C = 1, 

109.9  ft/sec 

M = 

0.0039565 

e = 

53°.  0 

True  Distance 

True  Time 

Time  Data 

Time  Error 

Mike 

Tj  (feet) 

r^c  (sec; 

t($ec) 

r^/c-t  (sec) 

1.) 

W6 

16,898.89 

15.226 

15.195 

.031 

2.) 

W7 

19,087.77 

17.198 

17.169 

.029 

3.) 

W8 

21,962.93 

19.788 

19.764 

.024 

True  (tj-t^)  = 1.972  sec 

(t  -tj  Data  = 1*974  sec  (long) 

True(t3-t^)  = 4.563  sec 

(t3-t^)Data=  4.569  sec  (long) 

True  $2  “ 2,188.88 

TrveS3=  5,064.04' 

Data  $2  = 2,190.94' | ^ 

Data  S3  = 5,071.13'  ) ® 


Wind  $2  = 
Wind  S3  = 


2,201.49') 
5,092.23  I v/ind 


16,786.16' 

16,697.83' 


Tro*  4.5-.*o 


0.6674  low 

1. 194  low 
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Day:  11/8/74(312) 
Even^:  92 
Gun:  2 

t = 17:20:00.178  MST 


B6  (31) 

B7  (32) 

B8  (33) 
d = 4,430.30' 


W = 3. 9 kfs  587  ft/sec) 
from  266  .4  at  25'elev. 

M = 0.0059368 
6 = 26°.  8 


T=  52. 1°F  (11.2°C)at26'elev. 
C = 1, 109. 5 ft/sec 


Mike 

True  Distance 
(feet) 

True  Time 
r^/c  (sec) 

Time  Data 
t(sec) 

Time  Error 
r^/c-t  (sec) 

1.) 

B6 

9,405.76 

8.477 

8.472 

+.005 

2.) 

B7 

12,915.59 

11.641 

11.645 

-.004 

3.) 

B8 

16,864.62 

15.200 

15.220 

-.020 

True  (t2-f]) 

= 3. 163  sec 

(»2”fi) 

= 3.173  sec  (long) 

True  (fj-t^) 

= 6.723  sec 

(t3-ti)  Data 

= 6.748  sec  (long) 

True  = 

3,509.83* 

True  $3  = 

7,458. 86* 

Data  S-  = 

3,520.44') 

2 

Data  Sg  = 

7,486.91')  o 

Wind  = 

3,543.92') 

2 

> R . j ~ 

Wind  = 

7,533.86')  wind 

8,954.88' 

8,536.40' 


:i3t3 

* w!  I 


4.7n  low 


9.  244  low 


r 


r 


149 


1.) 

2.) 

3.) 


Doy: 

11/8/74 

W6  (26) 

Event 

: 92 

W7  (34) 

Gun: 

2 

W8  (37) 

t = 
o 

17:20:00.178  MST 

d = 4,430.33' 

W = 

3.9  kts  (6.587  ft/sec) 

T = 52 

. 1°F  (11.2°C)  at  26* 

elev. 

from  266  .4  at  25'  elev. 

c = 1, 

1 09. 5 ft/ sec 

M = 

0.0059368 

6 = 

116°.  8 

True  Distance 

True  Time 

Time  Data 

Time  Error 

Mike 

r2  (feet) 

r^c  (sec) 

t(sec) 

r^/c-t  (sec) 

W6 

16,898.89 

15.231 

15.100 

.131 

W7 

19,087.77 

17.204 

17.073 

.131 

W8 

21,962.93 

19.795 

19.679 

.116 

True  = 1.973  sec 


Dq^o  = 1.973  sec  (no  error) 


True  = 4.564  sec 


(t3-f^)  Data  = 4.579  sec  (long) 


Tvoe  Jfj/o 

Vj-  n'c 
4-4!? 


True  $2  = 2,188.88' 


True  $2  - 5,064.04' 


Data  S-  = 2,189.04' 


Data  = 5,080.40' 


16,394.93'  ^ 2.98i  low 


Wind  $2  = 2,177.18' 


Wind  “ 5,056.68' 


wind 


16,492.38'  — 2.41^  low 


r 
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Doy: 

11/15/74  (319) 

A3  (14) 

Event: 

95 

A2  (13) 

Gun: 

1 

A1  (12) 

t » 

09:15:00.087  MST 

d ='4,430.70 

W = 2.9  kts  (4.898  ft/s«c  ) T “ 5l.4°F  (10.8°C)  at  26'  elov. 

from  2ll°.4  at  25'olev.  q „ 1, 108.7  ft/$oc 

M = 0.0044178 
0 = 261°.8 


Mike 

True  Distance 
(feet) 

True  Time 
r^/c  (see) 

Time  Data 
t($ec) 

Time  Error 
r^/c-\  (sec) 

A3 

9,410.51 

8.488 

8.456 

^032 

A2 

12,941.68 

11.673 

11.697 

-.024 

A1 

16,901.73 

15.245 

15. 280 

-.035 

True  “ 3.185  sec 

(t^-t^)  Data  * 3-241  sec  (long) 
True(t3-t,)  “ 6.757  sec 
(t  -t.)  Data  •»  ^-324  sec  (long) 

Tvet  = 

a»5*.o 

True  $2  “ 3,531. 17' 
True  Sg  * 7,491. 22' 
Data  $2  ® 3,593.30' 
Data  Sg  “ 7,565.77' 

j'o  - 

10,344.88' 

•*"  9. 93't'  high 

Wind  $2  “ 3,590.51' 
Wind  “ 7,560.  19'  i 

1 R - 

1 wind 

10,403.32' 

^ 10. 6<  high 

r 


Day: 

11/15/74  (319) 

B6  (31) 

Event 

: 95 

B7  (32) 

Gun: 

1 

B8  (33) 

t = 
o 

09:15:00.087  MST 

d = 4^430.30' 

W = 

2. 9 kts  (4. 898  ft/sec ) 

T = 51 

.4°F  (10.8°C)  at  26'elev 

from  211°.  4 at  25'  elev 

C = 1, 

108. 7 ft/sec 

. 

M = 

0.0044178 

0 = 

GO 

O * 

00 

• 

True  Distance 

True  Time 

Time  Data 

Time  Error 

Mike 

(feet) 

r /c  (sec) 

1 » 

t(sec) 

r^/c-t  (sec) 

1.) 

B6 

9,385.47 

•8.465 

8.478 

-.013 

2.) 

B7 

12,891.24 

11.627 

11.641 

- .014 

3.) 

B8 

16,838.64 

15.188 

15,210 

- .022 

TrueO^-t^)  = 3.162  sec 

Data  = 3.163  sec  (little  long) 

True  = 6.722  sec 

“ 6.732  sec  (long) 

True  $2  = 3,505.77' 

TrueSg  = 7,453.17' 

Data  $2  = 3,506.82'  1 ^ ^9  , 

Data  Sg  = 7,463.77'  ( ® 

Wind  = 3,509.61'  | =8 

Wind  S,  = 7,469.35'  ) 


Trot 

Y,-a.3*.o 


= 9,044.97'  3.63^  low 


R . .=  8,995.90' 
wind 


4. 15^  low 


[ 

Day;  11/15/74(319)  W6  (26) 

Event:  95  W7  (34) 

Gun:  1 W8(37) 

I t = 09:15:00.087  MST  d = 4,430.33' 


W = 

2.9  kts  (4.898  ft/sec) 

T = 51 

. 4°F  (10. 8°C)  at  26'  elev 

from  2l  1°.4  at  25'  elev 

C = 1, 

108.7  ft/sec 

M = 

0.0044178 

0 = 

00 

• 

o 

True  Distance 

True  Time 

Time  Data 

Time  Error 

Mike 

fj,  (feet) 

r^/c  (sec) 

t(sec) 

r^/c-t  (sec) 

1.) 

W6 

16,924.14 

15.265 

15.273 

-.008 

2.) 

W7 

19,109.90 

17.236 

17.253 

-.017 

3.) 

W8 

21,981.98 

19.827 

19.854 

-.027 

True  (t^-t^)  = 1*971  sec 

(t2~t^)  Data  = 1.980  jg^  (long) 

True  (»3-»^)  = 4.562  sec 

(t3-t^)  Data  = 4.581  sec  (long) 

True  $2  - 2,185.76' 

TrueSg  = 5,057.84' 

Data  $2  = 2,195.23'  1 ^ 

Data  Sg  = 5,078.95'  ) o 
Wind  $2  = 2,175.86'  j ^ 
Winds.  = 5,040.21'  I wind 


16,774.38' 

16,936.10' 


= ‘V^’3. 


T 
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Day:  11/18/74  (322) 
Event:  4 1 
Gun:  2 

06:34:00.293  MST 


A3  (14) 

A2  (13) 

A1  (12) 
d = 4,430.70* 


W = 1.7  kts  (|.  87  ft/sec)  f = 45. 8°F  (7. 7°C)  at  26*  elev 

from  173  .7  at  25*  elev  q = 1,102.7  ft/sec 

M = 0.0026038 
0 = 299°.  5 


True  Distance 

True  Time 

Time  Data 

Time  Error 

Mike 

Tj  (feet) 

r^/c  (sec) 

t(sec) 

r^/c-t  (sec) 

1.) 

A3 

9,390.03 

8.515 

8.358 

.157 

2.) 

A2 

12,917.24 

11.714 

11.620 

.094 

3.) 

A1 

16,875.70 

15.304 

15.161 

.143 

True  (t^-tj)  = 3*199  sec 

(t^-t^)  Data  = 3.262  sec  (long) 

True  (tg-t^)  = 6.788  sec 

Data  = 6.803  sec  (long) 

True  $2  = 3,527.21* 

True  $2  = 7,485.67* 

Data  =3,597.01*  ) 

Data  $2  =7,501.67*  | o 

Winds  =3,602.69*  j 

Winds  =7,513.03*  / wind' 


14,406.17* 

14,261.88* 


•*-  53.4'ii^  high 


51.9^  high 


r 
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Day; 

11/18/74  (322) 

A6  (17) 

d = 
2 

1 ;5720.99' 

Event 

: 41 

A2  (13) 

“s' 

4,430.45' 

Gun: 

2 

A1  (12) 

D = 

2 2J51.44' 

t = 
o 

06:34:00. 293  MST 

■ 

W = 

1.7  kts  (2. 87  ft/sec) 

T = 45 

.8°F(7.7°C)  « 

at  26'  elev 

from  173°.  7 at  25'elev 

C = 1, 

102.7  ft/sec 

M = 

0.0026038 

0 = 

299°. 5 

True  Distance 

True  Time 

Time  Data  Time  Error 

Mike 

(feat) 

r^/c  (sec) 

t(sec) 

r^/c-t  (sec) 

1.) 

A6 

9,364.19 

• 8.492 

8.382 

.110 

2.) 

A2 

12,917.24 

11.714 

11.620 

.094 

3.) 

A1 

16,875.70 

15.304 

15.161 

.143 

9,624.28' 


tk 
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Day: 

11/18/74  (322) 

85  (30) 

d2  = ^l.ll’ 

Event 

: 41 

87  (32) 

d^=  4,430.36' 

Gun: 

2 

88  (33) 

D = d2  + d3=l 

^1.47' 

t = 
o 

06:34:00. 293  MST 

W = 

1.7  kts  (2.87  ft/sec) 

T = 45. 

,8°F  (7.7°C)  at  26'  elev 

from  173°. 7 at  25'  elev 

C = 1, 

102.7  ft/sec 

M = 

0.0026038 

t 

0 = 

119°  5 

True  Distance 

True  Time 

Time  Data 

Time  Error 

Mike 

(feet) 

r^/c  (sec) 

t(sec) 

ry'c-t  (sec) 

1.) 

85 

7,026.36 

6.372 

6.368 

.004 

2.) 

87 

12,915.59 

11.713 

11.693 

.020 

3.) 

88 

16,864.62 

15.294 

15.279 

.015 

True  = 5.341  sec 

Data  = 5.325  sec  (short) 

True  “ 8.922  sec 

(t3-t^)  Data  =8.911  sec  (short) 

True  $2  ~ 5/889.23 

True  Sg  = 9,838.26 

Data  S “ 5,871.88  | 

2 ) R 

Data  Sg  = 9,826.16'  | o 

Wind  = 5,860.52'  ) 

2 > R 

Winds  =9,801.12'  j wind 

3 


6,897.74* 


1.83^  low 


6,963.86' 


0.605^  low 
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Day: 

11/18/74 

B6  (31; 

Event 

: 41 

B7  (32) 

Gun; 

2 

B8  (33) 

t = 

0 

06:34:00. 293  MST 

d = 4,430.30' 

W = 

1.7  kts  (2.87  ft/sec) 

T = 45.8°F  (7 

'.7°C)  at  26'  elev 

from  173^.7  at  25'  elev 

C = 1,102.7 

ft/sec 

M = 

0.0026038 

e = 

119°  5 

True  Distance 

True  Time 

Time  Data 

Time  Error 

Mike 

(feet) 

T^c  (sec) 

t(sec) 

r^c-t  (sec) 

1.) 

B6 

9,405.76 

8.530 

8.517 

.013 

2.) 

B7 

12,915.59 

11.713 

11.693 

.020 

3.) 

B8 

16,864.62 

15.294 

15.279 

.015 

True  ~ 3.183  sec 

(^2-^^)  Data  = 3.176  (short) 

True  (f3-t|)  = 6.764  sec 

(t^-t^)  Data  = 6.762  (short) 

True  $2  = 3,509.83' 

True  $2  = 7,458.86' 

Data  $2  =3,502.18'  ) ^ 

Data  Sg  =7,456.46'  ( <> 

Wind  =3,496.50'  ) 

2 \ R 

Wind  S„  =7,445.10'  ) w 


'Troc  7,1.2 


9,054.25'  — 3.74<low 


9,153.54'  — 2.68ilow 
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Day: 

11/18/74  (322) 

W6  (26) 

Event 

; 41 

W7  (34) 

Gun: 

2 

W8  (37) 

t = 

0 

06:34:00. 293  MST 

d = 4,430.33' 

W = 

1.7  kts  (2. 87  ft/sec) 

T = 45. 

.8°F  (7.7°C)  at  26'  elev 

from  173° 7 at  25'  elev 

C = 1, 

102.7  ft/sec 

M = 

0.0026038 

e = 

209°  5 

True  Distance 

True  Time 

Time  Data 

Time  Error 

Mike 

(feet) 

r^c  (sec) 

t(sec) 

r^c-t  (sec) 

1.) 

W6 

16,898.89 

15.325 

15.260 

.065 

2.) 

W7 

19,087.77 

17.310 

17.247 

.063 

3.) 

W8 

21,962.93 

19.917 

19.864 

.053 

True  (12“^])  “ 1.985  sec 

Data  = 1.987  sec  (long) 
True  = 4. 592  sec 

(t^-ti)  Data  = 4.604  sec  (long) 
True  $2  = 2,188.88' 

TrueSg  = 5,064.04' 

Datb  $2  = 2,191.06'  ) ^ 

Data  $2  = 5,076.83'  ( ® 

^ Wind  $2  = 2,181.02'  j ^ 

Winds  =5,056.75'  ) wind 

w 


Tru.  4-s'.*o 

^3  = W 

16,613.59'  - 1.69^  low 


t 


I 


16,696.62 


1. 20^  low 
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T 


1 


1.) 

2.) 

3.) 


Day;  11/19^4  (323) 

Event:  13 

Gun:  1 

t = 06:36:59.968  MST 
o 

A3  (14) 

A2  (13) 

A1  (12) 
d = 4,430.70* 

W = 

5.9kts(9.97  ft/sec) 

T = 56.4°F  (13.5°C)at  26' 

elev 

from  274°  8 at  25'  elev 

C = 1,114.0  ft/sec 

M = 

0.008945 

e = 

198% 

True  Distance 

True  Time 

Time  Data 

Time  Error 

Mike 

Tj  (feet) 

r^/c  (sec) 

t(sec) 

r^/c-t  (sec) 

* 

A3 

9,410.51 

8.447 

8.325 

.122 

A2 

12,941.68 

11.617 

11.504* 

.113 

A1 

16,901.73 

15.172 

15.099 

.073 

ti 

i 

! 


True  (h  -f.)  = 3. 170  sec 
2 1 

(yt,)  Data  =3. 179  sec  (long) 
True(ytp  = 6.725  sec 
(t^-t,)  Data  = 6.774  sec  (long) 
True  $2  = 3/531.17' 

True  Sg  = 7,491.22' 

Data  =3,541.41'  1 
Data  = 7,546.24'  j *^0 
Wind  S = 3,503.80  | 

Wind  Sj  =7,471.03'  J "wind 


7,983.73' 

8,630.37' 


( * Traces  available  for  these  mikes) 


"Tvo*  0 

a\'.t 


15.2;^  low 
8.  29^  low 


I 


i 


"i  ri*r< 
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Day;  11/19/74  (323) 
Event:  13 
Gun:  1 

t = 06:36:59. 968  MST 
o 


A6  (17) 

A7  (18) 

A8  (19) 
d = 4,430.39' 


W = 5.9  kfs  (9.97  ft/sec) 
from  274°  8 at  25'  elev 


T=  56.4°F(13.5°C)at26'  elev 
C = 1, 114.0  ft/sec 


M = 0.008945 
e = 341°6 


True  Distance 

True  Time 

Time  Data 

Time  Error 

Mike 

r^  (feet) 

r^/c  (sec) 

t(sec) 

r^  /c-t  (sec) 

1.) 

A6 

9,345.14 

8.389 

8.079* 

.310 

2.) 

A7 

12,862.16 

11.546 

11.116* 

.430 

3.) 

AS 

16,816.25 

15.095 

14.529 

.566 

True  (t^-t^  =3. 157  sec 

^^2*^1^  Data  = 3.037  sec  (short) 

True  (»3-t|)  = 6.707  sec 

(to-t^)  Data  = 6.450  sec  (short) 

True  $2  = 3,517.02' 

True  $3  = 7,471.11* 

Data  $2  = 3,383.22'  ) ^ 

Data  = 7,185.30'  ( o 

Wind  =3,420.82'  | 

2 > R 

Winds.  =7,260.51’  / wind 


12,558.53' 

11,872.09' 


7,1.  H- 


-►  34.4^  high 
^ 27.0^  high 


r 


( * Traces  available  for  these  mikes) 
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Day:  11/19/74  (323) 

Evenf:  13 

Gun:  1 

t = 06:36:59.968  MST 
o 


B3  (28) 

82  (27) 

81  (26) 
d = 4,430. 13' 


\V  = 5.9  kts  (9.97  ft/sec) 
from  274°  8 at  25'  elev 

M = 0.008945 
0 = 161°6 


j = 56.4°F  (13.5°C)  at  26'  elev 
C = 1, 1 14.0  ft/sec 


Mike 

True  Distance 
r^  (feet) 

True  Time 
(sec) 

Time  Data 
t(sec) 

Time  En 
r^^/c-t 

1.) 

83 

9,451.63 

8.484 

8.545 

-.061 

2.) 

82 

12,970.57 

11.643 

11.643* 

.000 

3.) 

81 

16,924.14 

15.192 

15.152 

+.040 

Troe(t2-t^)  = 3. 159  sec 
^^2"^1^  Data  = 3.098  sec  (short) 
True(t3-t^)  = 6.708  sec 
(t^-t^)  Data  = 6.607  sec  (short) 
True  $2  = 3,518.94' 

True  Sg  = 7,472.51 
Data  ~ 3,451.17'  \ 

Data  Sg  = 7,360.20'  | o 
Wind  $2  ~ 3,413.57  | 

Winds..  = 7,285.00'  | *^wind 


Trvt  :\3!x 
3 


9,719.90' 

10,358.84' 


9.6D?  high 


( * Trace  available  for  this  mike) 
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i 

% 

i Day:  {353) 

Event:  \ 3 

T Gun:  \ 

MST 


\:)4-C:iO 

AG  cn) 

£2.  Ml) 

c)  = 


I 


i 

i 


W = S.'lKts  (.^.<11  ft(sec)  T = (13.  5"'^:)  2t' tV\/. 

c = 1,114".  0 ?v/s<c. 

M = O.oo?«1tfy 

e = tx 


Mike 

True  Distance 

Tj^  (feet) 

True  Time 
r^/c  (sec) 

Time  Data 
t(sec) 

Time  Error 
rj/c-t  (sec) 

M- 

3,\As-.?1 

A-m* 

A G 

A,34s-.l'f 

S.oil* 

.■310 

1%  00  1 . a? 

n.osi 

/t.  ‘/V? 

» 6 fO 

TrueCtj-t^)  =?.  Sts' set 
^^2**1^  Data  = S.  3®^  sec  C$V»>^) 
True  (^2"!'^)  “ ^33  sec 

(ts-ti)  Data  = 13.((>n:istt  (sWi) 
True  $2  = :^s' 

True  = lSr,?5S.3l' 

Data  $2  = S.^OT-tt’  ) ^ 

Data  Sg  = \S^:i3£>.tl'  ( ° 

Wind  $2  = S', ‘I It. 31 ' ) ^ 

WindSg  = lS,‘f£)t.03'/ 


H-,A44.oi' 

4,511.13' 


True  1/^  = 1^0 

Y,  - h:3 

I3*  1 


" 4^-g;  wv 


( Trrfcs  <vvwU\>\e  W Wvesr  n>iH«-s) 


i 
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Event:  \3 
Gun:  \ 


T)5U^)  \a,^531.M?' 

AT  Cl?)  ^3= 

EX  (11) 


W =5.‘lK^  I'l-'n  ^/yt)  T “ 5(..q‘f  AVat'tUv. 


i 

M = 

e = 

fvov>^  n!f.7  J^tVV. 

0.  ooT'iHS' 

= i> 

1 L 4-.  0 -ft  /;-cc . 

True  Distance 

True  Time 

Time  Data 

Time  Error 

Mike 

r^  (feet) 

r^/c  (sec) 

t(sec) 

r^/c-t  (sec) 

.)  ^sr 

3,uo.  ax 

— . oo3 

.)  i^T 

\x,ux.i(, 

u.  lie* 

+ . T30 

.)  El 

l<\,ool.  a? 

n.  65-7 

It. 41^7 

+ . t/o 

True  (t^-t^)  = ^.THS'sec 

Data  = ?.i\2-9ec  CsV»*-'t) 

True  ~ \4-.aS’t  kc 

(ta-t^)  Data  = I'i.  (^43$et  CiV»'^t) 

True  $2  = 

True  $2  = \S,it\.f>{,' 

DataSj  = ^.a5l.57'  ) ^ 
DatoSg  = \S,m.3t>')  o 
Wind  $2=  ‘\,35<).?o'l  ^ 

WindSa  = 1S-,3'^8.4V7  wind 


S,?U.57' 


Trot  ^ ^.“s- 

% - i^rc 

3 


- 

66.7  % Uia.W 


(*TrKtJ  (SVM\*\)\t  Wtit  vw'Vt;.^ 


Day; 

Event:  l3 
Gun;  1 

t = W\ST 


W = 5,1  Kts  (q.11 

^'fcw\  iT4,  ? (xt  is’fcVv. 

M = O.00?rf5' 

e = 


^Ic  (33) 

El  (30 
E5  M) 

4,430. 3i' 


T = 5t,.4^F  (13.5‘C)  <aV  ;it'e,Vv. 

C = 1,1  H.O  tt/$«c 


Mike 

True  Distance 
r^  (feet) 

True  Time 
r^/c  (sec) 

Time  Data 
t($ec) 

Time  Error 
r^/c-t  (sec) 

1.) 

Et 

lk,n«.(rf 

is-.  n5 

14. 134 

. 1% 

2.) 

ET 

14,03'f.  14 

n-  on 

I6.n?4 

3.) 

E« 

14.  673 

14-461 

True  (^2"^]^  “ 

(t2-t^)  Data  = ^.OSO  stc- 
True  (^3“^|)  “ H.SS?  s«- 
(t^-t^)  Data  = 4.1  S<c  (loa) 
True  $2  = 

True  $2  = 5,on.SS' 

DataS2  = 3,a?3.V  I ^ 

Data  Sg  = S-,H0.3\'  I o 

WlndS2=  A,311.\V  j ^ 

Wind  Sg  = ^,  335. »')  wind 


t^uc  zn 

r,  - 


R„  = \'5,%K.c^^'  - 

' .„H=  ♦ 


I2  Jeiv 

S.  03  ^ l»iv 
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Day;  U/}9/74  {323) 
Event;  34 
Gun:  1 

t = 07:34:00.389  MST 
o 


\V  = 5.4  kts  (9. 12  ft/sec) 
from  137^23  at  25’  elev 

M = 0.008180 
e = 335° 9 


A3  (14) 

A2  (13) 

A1  (12) 
d = 4,430.70' 


T = 57.2°F  (14.0°C)  at  26’  elev 


C = 1,115.0  ft/$ec 


Mike 

True  Distance 
(feet) 

True  Time 
r^/c  (sec) 

Time  Data 
t(sec) 

Time  Error 
r^/c-t  (sec) 

1.) 

A3 

9,410.51 

8.440 

8.445* 

-.005 

2.) 

A2 

12,941.68 

11.607 

11.698* 

-.091 

3.) 

A1 

16,901.73 

15.159 

15.329 

-.170 

True  (t^-t^)  = 3. 167  sec 

(tj-t^)  Data  = 3.253  sec  (long) 

True  (tg-tj)  = 6.719  sec 

(t^-t^)  Data  = 6.884  sec  (long) 

True  $2  = 3,531. 17’ 

True  $2  = 7,491.22’ 

Data  $2  ~ 3,627.10’  ) 

Data  Sg  = 7,675.66'  | o 

Wind  S = 3,660. 18’  | 

2 > R " 

Winds..  = 7,741.83’  ) wind 


7,898.92’ 

7,260.43' 


T-cot  1(3.0 


16. 1^  low 
22.84  low 


( *Traces  avaiicble  for  these  mikes) 


I 


Day:  11/19/74  (323) 
Event:  34 
Gun:  1 


t = 07:34:00.389  MST 
o 


A6  (17) 

A7  (18) 

A8  (19) 
d = 4,430.39' 


W = 

5.4  kts  (9. 12  ft/sec) 
from  13^3  at  25'  elev 

T = 57 

C = ^ 

.2°F  (14.0®C)  at  26' 

1 15.0  ft/sec 

elev 

M = 

0.008180 

e = 

204°  1 

Mike 

True  Distance 
r^  (feet) 

True  Time 
r^/c  (sec) 

Time  Data 
t(sec) 

Time  Error 
r^/c-t  (sec) 

1.) 

A6 

9,345. 14 

8.381 

8.238* 

.143 

2.) 

A7 

12,862.16 

11.536 

11.384* 

.152 

3.) 

A8 

16,816.25 

15.082 

14.806 

.276 

True  (t^-t^)  = 3. 154  sec 


Data  = 3.146  sec  (short) 


True  ~ 6.701  sec 


(to-t^)  Data  = 6.568  sec  (short 


Tfv.c 


True  $2  =3,517.02' 


True  $2  =7,471.11' 


Data  $2  =3,507.79' 
Data  S =7,323.32' 


Wind  $2  = 


Wind  S, 


=3,474.71'  ) 
=7,257.16'  ) 


16,629.10' 

17,445.83' 


^ 77. 9i  high 


( * Traces  available  for  these  mikes) 


— ^ I I , , II  ...  - v«.».  ■ . - ■ ^ — .1.  ■ ^_., 
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! 


' Day:  11/19/74  (323) 

85  (30) 

d2  = 8,860.58* 

! Event:  34 

83  (28) 

d2  = 4,429.34* 

Gun:  1 

t 

82  (27) 

D=  13,289.92* 

t = 07:34:00.389  MST 
o 

mn\ III  iij  ■ 


( 


\V  =5.4  lets  (9. 12  fl/sec) 
from  13^.3  at  25'  elev 


M = 0.008180 
0 = 24°.! 


T=  57.2°F  (14.0°C)at26' elev 
C = 1,115.0  ft/sec 


True  Dislance  True  Time  Time  Data  Time  Error 


Mike 

r^^  (feet) 

r^/c  (sec) 

t(sec) 

r,  /c-t  (sec) 
1. 

1.)  B5 

7,016.80 

6.293 

6.311 

-.018 

2.)  B3 

9,451.63 

8.477 

8.591 

-.114 

3.)  S2 

12,970.57 

11.633 

11.793* 

->.160 

True  ~ 2. 184  sec 

(t2-ti)  Data  = 2.280  sec  (long) 
True  (tg-t^)  = 5. 340  sec 
(t^-t^)  Data  = 5.482  sec  (long) 
True  $2  = 2,434.83* 

' TrueSg  =5,953.77* 

Data  =2,542.20*  ) 

I 2 ) R ■■ 

DataSg  =6,112.43  | o 

Wind  =2,608.36*  1 

I Wind  $2  =6,211.67*  f wind 

i 


31.0 
30*0 
= ^<3 

6,783.77*  — 3.32^  low 

6,628.94*  ^ 5.531b  low 


( * Trace  available  for  this  mike) 
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Day:  11/19/74  (323)  D4  (21) 

Evenf:  34  A6  (17) 

Gun:  1 E2(11) 


t = 07:34:00.389  MS T d = 9,906.40' 


w = 

5.4  kts  (9. 12  ft/sec) 
from  137°.  3 at  25'  elev 

T = 

C = 

57.2°F  (14.0°C)  at  26’  elev 

1, 115.0  ft/sec 

M = 

0.008180 

0 = 

o 

o 

• 

Mike 

True  Distance 

Fj  (feet) 

True  Time 

r /c  (sec) 
1 

Time  Data 
t(sec) 

Time  Error 
r^/c-t  (sec) 

D4 

3,145.89 

2.821 

2.815* 

.006 

A6 

9,345.14 

8.381 

8.238* 

.143 

E2 

19,001.28 

17.042 

16.806 

.236 

True  (t^-t^)  = 5. 560  sec 
^^2"h^  Data  = 5.423  sec  (short) 
True  (tg-t^)  = 14. 220  sec 

(to-t^)  Data  = 13.991  sec 

True  $2  =6,199.25’ 

True  $2  = 1 ^55.39' 

Data  $2  =6,046.65'  ) ^ 

Data  =15,599.97’  ) <» 

Wind  =5,995.32'  ) 

2 5 R . 

Winds  =15,497.32')  wind 


~Wy«-  ^ = ^^6 
%=  v>:  s 
u'.t 

3,712.70’  ^ 18.0^  high 
3,991.57'  ^ 26. 9^  high 


( * Traces  available  for  these  mikes) 


— f 
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Day;  11/19/74  (323) 
Event;  34 
Gun:  1 

t = 07:34:00.389  MST 


E6  (33) 
E7(36) 

E8  (44) 

d = 4,430.32' 


\V  = 5.4  lets  (9. 12  ft/sec) 
from  13^.3  at  25'  elev 


M = 0.008180 
0 = 114®.! 


T = 57.2°F  (14.0®C)  at  26'  elev 
C = 1, 1 15.0  ft/sec 


Mike 

True  Distance 
(feet) 

True  Time 
r^/c  (sec) 

Time  Data 
t(sec) 

Time  Error 
r^^'c-t  (sec) 

1.) 

E6 

16,838.64 

15.102 

14.897 

.205 

2.) 

E7 

19,034.19 

17.071 

16.864 

.207 

3.) 

E8 

21,916.19 

19.656 

19.555 

.101 

True  (t^-tp  = 1.969  sec 

(t^-t^)  Data  =1.967  sec  (short) 

True(t3-t^)  =4.554  sec 

(tg-t^)  Data  =4.658  sec  (long) 

True  $2  = 2,195.55' 

True  Sg  =5,077.55' 

Data  S.,  =2,193.21'  | 

Data  Sg  =5,193.67'  | o 

Wind  =2,178.41' 

2 

Wind  S = 5,164 


,4V  . 

.or  I 


wind 


13,565.56' 

13,675.31' 


True 


— 19.451  low 


^ 18.85g  tow 
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Day:  12/Z^4  (336) 
Event:  9 
Gun:  2 

t = 06:32:02.083  MST 


A3  (14) 

A2  (13) 

A1  (12) 
d = 4,430.70' 


W=  1.8  kts  (3.04  ft/sec)  T=  36.3°F  (2.4°C)  at  26'  elev 

from  200°.  1 at  25'  elev  ^ ^ ],092. 2 ft/sec 

M = 0.0027835 
0 = 273°  1 


Mike 

True  Distance 
(feet) 

True  Time 
r^c  (sec) 

Time  Data 
t(sec) 

Time  Error 
r^/c-t  (sec) 

A3 

9,390.03 

8.597 

8.552 

+.  045 

A2 

12,917.24 

11.827 

11.882 

-.055 

A1 

16,875.70 

15.451 

15.507 

-.056 

True(t2-t^)  = 3. 229  sec 

(t^-ti)  Data  = 3.330  sec  (long) 

True  (tg-tp  = 6. 854  sec 

(t_-t.)  Data  = 6.955  sec  (long) 

T'cvt  53.  \ 

4iwiva  1 

True  $2  = 3,527.21' 
True  $2  = 7,485.67' 
Data  $2  = 3,637.03'  | 
Data  = 7,596.25'  I 

R 

1 ° 

12,438.57' 

32. 5^  high 

Wind  $2  = 3,637.70'  ] 
Wind  Sg  = 7,597.58'  ) 

! R . .= 

1 wind 

12,422.72' 

— 32. 3<  high 
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Day: 

12/2/74  (336) 

86  (31) 

Event 

: 9 

87  (32) 

Gun: 

2 

88  (33) 

t = 

0 

06:32:02.083  MST 

d = 4,430.30' 

W = 

1.8  kts  (3.04  ft/sec) 

T = 36.3°F  (2.4°C)  at  26'  elev 

from  200°.  1 at  25'  elev 

C = 1,092.2 

ft/sec 

M = 

0.0027835 

0 = 

93°.  1 

True  Distance 

True  Time 

Time  Data 

Time  Error 

Mike 

r^  (feet) 

r^/c  (sec) 

t(sec) 

r^/c-t  (sec) 

1.) 

B6 

9,405.76 

8.612 

8.519 

.093 

2.) 

87 

12,915.59 

11.825 

11.703 

.122 

3.) 

88 

16,864.62 

15.441 

15.311 

.130 

True  (*^2"^])  ~ 3.214  sec 

Dafa  = 3.184  sec  (short) 

True  ~ 6.829  sec 

(t^-t^)  Data  = 6.792  sec  (short) 

True  $2  = 3,509.83' 

True  $2  =7,458.86' 

Data  S.,  = 3,477.56'  | 

DataSg  =7,418.22'  ) <» 

Wind  =3,476.89'  ) 

2 } R . j 

Wind  S = 7,416.89'  j wind 

o 


9,082.24' 


True  v(>j=  ;^3.*3 

M?>=  '3.3' •i- 

3.44^  low 


f 

4 


9,093.29' 


3.32^  low 


172 


Day: 

12/2/74  (336) 

E6  (33) 

Event: 

9 

E7  (36) 

Gun: 

2 

E8  (44) 

t = 
o 

06:32:02.083  MST 

d = 4,430.32' 

W = 

1.8  kts  (3.04  ft/sec) 

T = 36 

.3°F  (2.4°C)  at  26' 

elev 

from  200°.  1 at  25'  elev 

C = 1/ 

092. 2 ft/sec 

M = 

0. 0027835 

e = 

176°.  9 

True  Distance 

True  Time 

Time  Data 

Time  Error 

Mike 

r^  (feet) 

r^c  (sec) 

t(sec) 

r^/c-t  (sec) 

E6 

16,864.62 

15.441 

15.311 

.130 

E7 

19,057.40 

17.449 

17.385 

.064 

E8 

21,936.54 

20.085 

19.985 

.100 

True  ~ 2.008  sec 

(f2*'lj)  Data  = 2.074  sec  (long) 

True  “ 4.644  sec 

(t^-t^)  Data  = 4.674  sec  (long) 

True  $2  = 2,192.78' 

True  = 5,071.92' 

Data  $2  = 2,265.22'  ) ^ 

Data  $2  = 5,104.94'  | <> 

Wind  = 2,252.91'  ) 

2 > R 

Winds.  = 5,080.31'  ) wind 

w 


44*1 


20,415.58' 

20,537.46' 


21. 8^  high 
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Day;  }2/2/74  (336) 

A3  (14) 

Event:  94 

A2  (13) 

Gun:  2. 

A1  (12) 

t = 09:44:00.075  MST 
o 

d = 4, 430, 70 

W = 1.4  kts  (2.365  ft/sec) 
from  281°.  7 at  25'  elev 

M = 0.002152 
e = 191°.  5 


T = 42.3°F  (5.7°C)  at  26'  elev 
C = 1,098.7  ft/sec 


Mike 

True  Distance 
r2  (feet) 

True  Time 
r^c  (sec) 

Time  Data 
t(sec) 

Time  Error 
r^/c-t  (sec) 

1.) 

A3 

9,390.03 

8.546 

8.476 

+.070 

2.) 

A2 

12,917.24 

11.757 

11.751 

+.006 

3.) 

A1 

16,875.70 

15.360 

15.378 

-.018 

True(t2-t^)  = 3.210  sec 

~ 3. 275  sec  (long) 
True  (tg-t^)  = 6.813  sec 
(*2“*'j)  Dafa  = 6.902  sec  (long) 
True  $2  = 3,527.21' 

True  $2  = 7,485.67' 

Dota  $2  =3,598.24'  ) ^ 

Data  =7,583.23'  ( ° 

Wind  =3,588.90'  j ^ 

Wind  =7,564.54'  ) 


9,892.04' 

10,084.74' 


T'fv.e  ‘XiA 

Hi= 


5.35%  high 
7.40^  high 


Shock  Wave  Influence 
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Day:  12/2/74  (336) 
Event:  94 
Gun:  2 

t = 09:44:00.075  IAST 
o 


A6  (17) 

A7  (18) 

A8  (19) 
d = 4,430.41' 


! 

! W = 1.4  lets  (2.365  ft/sec) 

from  281°.  7 at  25'  elev 

M = 0.002152 
e = 348°.  5 


T = 42.3°F  (5.7°C)  at  26'  elev 
C = 1,098.7  ft/sec 


Mike 

True  Distance 
r^  (feet) 

True  Time 
r^/c  (sec) 

Time  Data 
t(sec) 

Time  Error 
r^/c-t  (sec) 

1.)  A6 

9,364. 19 

8.523 

8.462 

.061 

2.)  A7 

12,885.60 

11.728 

11.694 

.034 

3.)  A8 

16,841.52 

15.329 

15.295 

.034 

True  “ 3.205  sec 

(t^-t^)  Data  = 3.232  sec  (long) 
True  “ 6.806  sec 

(ts-ti)  Data  = 6.833  sec  (long) 


True  $2  = 3,521.41' 
True  Sg  = 7,477.33' 
Data  $2  = 3,551.00' 
DafaSg  = 7,507.42' 

1 

= 10,008.03' 

6.88^  high 

Wind  = 3,560.34'  | 

! R . 

9,825.09' 

^ 4.92^high 

Wind  S = 7,526.11'  j 

( wind 

'V,=  A 3°.  I 


Shock  Wave  Influences 


M = 0.002152 


True  Dtstunce 


Wind  S 
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Day;  12/2/74  (336) 
Event;  94 
Gun;  2 

t = 09;44:00.075  MST 
o 


E3  (19) 

E2(ll) 

El  (8) 

d = 4,430.33' 


W = 

1.4  kts  (2.365  ft/sec) 
from  281°. 7 at  25'  elev 

T = 42.3°F  (5.7°C)at  26' 
C = 1,098.7  ft/sec 

elev 

M = 

0.002152 

e = 

281°. 5 

i 

1 

i 

1 

1 

Mike 

True  Distance 
rj  (feet) 

True  Time 
r^c  (sec) 

Time  Data 
t(sec) 

Time  Error 
rVc-t  (sec) 

' 1.) 

E3 

16,841.52 

15.329 

15.295 

.034 

2.) 

E2 

19,023.42 

17.314 

17.271 

.043 

3.) 

El 

21,895.15 

19.928 

19.780 

. 148 

True  (tj-t^)  = 1.986  sec 

(f2-t^)  Data  = 1 . 976  sec  (short) 

True  ~ 4.600  sec 

(t^-t^)  Data  = 4.485  sec  (short) 

True  $2  = 2,181.90' 

TrueSg  = 5,053.63* 

Data  $2  “ 2,171.03'  ) 

Data  Sg  = 4,927.67'  ) «> 

Wind  = 2,172.93'  1 ^ 

2 ) R ~ 

Winds.  = 4,931.47'  ) wind 


20,833.37' 

20,815.48' 


T>rvt  Vs-.O 

H7'C 
- 47. V 


23.74  high 


23. 6<  high 


Shock  Wove  Influence 


S E 2 

CO  tN  •— 
UJ  UJ  Ui 


>» 

p* 

p* 


«/> 

«o 

Ns 

o 

• 

o 
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Day:  12/2/74  (336) 
Event:  94 
Gun:  2 

t = 09:44:00.075  MST 
o 


W6  (26) 

W7  (34) 

W8  (37) 
d = 4,430.33' 


W = 1.4  kts  (2.365  ft/sec) 

T = 42.3°F  (5.7°C)at  26' 

elev 

from  281°.  7 at  25'  elev 

C = 1 

,098.7  ft/sec 

M = 0.002152 

0 = 101°.5 

True  Distance 

True  Time 

Time  Data 

Time  Error 

^^Hce  r^  (feet) 

r^/c  (sec) 

t(sec) 

r^/c-t  (sec) 

1.) 

W6  16,898.89 

15.381 

15.389 

00 

o 

o 

1* 

2.) 

W7  19,087.77 

17.373 

17.376 

-.003 

3.) 

W8  21,962.93 

19.990 

19.979 

+.011 

True  (t^-t^)  = 1.992  sec 

' 

^s.’o 

^^2*^1^  Data  = 1.987  sec  (short) 

iruc  '*2 

(t3-ti)  Data  = 4.590  sec  (short) 
True  $2  = 2,188.88' 

True  $2  = 5,064.04' 

Data  $2  = 2,183.12'  ) 

Data  = 5,043.03'  | ® 

Wind  $2  “ 


- HS‘.3 


Wind  S, 


= 2,181.22'  ) 
= 5,039.23'  ) 


wind 


17,254.63' 

17,270.68' 


2.1H  high 
2. 20^  high 


r 
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Day:  12/Z/74  (336) 
Event:  176 
Gun:  2 

t = 13:44:00. 180  MST 


B6  (31) 

B7  (32) 

B8  (33) 
d = 4,430.30' 


W = 2.4  kts  (4.05  ft/sec)  T = 54. 1°F  (12. 3 C)  at  26' elev 

from  268°.3  at  25'elev  ^ ^ 

M = 0.0036465 
e = 24®.  9 


True  Distance 

True  Time 

Time  Data 

Time  Error 

Mike 

(feet) 

r^c  (sec) 

t(sec) 

r^c-f  (sec) 

1.) 

B6 

9,405.76 

8.461 

8.452 

.009' 

2.) 

B7 

12,915.59 

11.619 

11.611 

.008 

3.) 

B8 

16,864.62 

15.171 

15.158 

.013 

True  (t^-t^)  = 3. 157  sec 
(t^-t^)  Data  = 3. 159  sec  (long) 
True  (tg-tp  = 6.710  sec 
(t.,-t,)  Data  = 6.706  sec  (short) 

w I 


True  $2  = 3,509.83' 

True  $2  = 7,458.86' 

Data  $2  =3,511.54'  1 

Data  Sg  = 7,454.39'  j 

Rq  = 9,678.67' 

-►  2.90^  high 

Wind  $2  = 3,526.19'  1 

D 9,409.93' 

0.0443^f  high 

Wind  S ~ 7,483.70*  j 

w 

wind  (4.  17'  off) 

True. 
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I 


Doy;  1^2/74  (336) 

W6  (26) 

Event:  176 

W7(34) 

Gun:  2 

W8  (37) 

t = 13:44:00. 180  MST 
o 

d = 4, 430. 33 

I 


W = 2. 4 Icfs  (4.05  ft/sec)  T = 54. 1°F  (12. 3®C)  at  26'  elev 

from  2680.3  at  25' elev  ^ ^ 1, 111.6  ft/sec 

M = 0.0036465 
e = 114°.9 


Mike 

True  Distance 
r^  (feet) 

True  Time 
r^c  (sec) 

Time  Data 
t(sec) 

Time  Error 
r^/c-t  (sec) 

1.) 

W6 

16,898.89 

15.202 

15.186 

.016 

2.) 

W7 

19,087.77 

17.171 

17.168 

.003 

3.) 

W8 

21,962.93 

19.758 

19.731 

.027 

True  (tj-t^)  = 1.969  sec 
^^2~*1^  Data  = 1.982  sec  (long) 
True  (tg-t^)  = 4.556  sec 
(t3-t^)  Data  = 4.545  sec  (short) 
True  $2  = 2,188.88' 

True  $2  = 5,064.04' 

Data  S = 2,203.19'  | 

Data  Sg  = 5,052.22'  ( <> 

WindS^  = 2,196.39'  1 

2 5 R . j 

Wind  S = 5,038.62'  j wind 


'TVut  4-5- *0 


18,145.99'  — 7.38^  high 
18,205.91'  ^ 7. 73^  high 


; 


I : 

' Day:  12/3/74  (337?  A3  (14)  • 

Event:  91  A2  (13) 

Gun:  2 

d = 4,430.70' 

‘ 10:13:59.771  MST 


“3.8  kts  (6.42  ft/sec)  T “ 47.8°F  (8.8°C)  at  26'  elev 

from  183°.  6 at  25'  elev  C = 1, 104.8  ft/sec 

M = 0.0058092 
e = 289°. 6 


I 


Mike 

True  Distance 
(feet) 

True  Time 
r^c  (sec) 

Time  Data 
t(sec) 

Time  Error 
T^c-t  (sec) 

1.) 

A3 

9,390.03 

8.499 

8.447 

+.052 

i 2.) 

t 

A2 

12,917.24 

11.692 

11.689 

+.003 

' 3.) 

A1 

16,875.70 

15.275 

15.349 

-.074 

True  (f2"^])  = 3. 193  sec 
(t2“t^)  Data  = 3.242  sec  (long) 
True  (^3"^^)  “ 6.776  sec 
(t3-t^)  Data  = 6.902  sec  (long) 
True  $2  = 3,527.21' 

True  Sg  = 7,485.67' 

Dahi  $2  = 3,581.76'  ) ^ 

Data  $2  = 7,625.33'  | ° 

Wind  $2  = 3,590.39'  ) 

Wind  S-  = 7,642.60'  ) 

w 


True  ^3 . 1 

= ^0. 1 

334.79'  ^ 21.9^  low 

183.18'  ^ 23.5^  low 


i 


k 


FT 


f 


'^Cv. 


r 
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, ! 


Day:  ^ 2/3/74  (337) 
Event:  91 
Gun:  2 

t = 10:13:59.771  MST 


i 

B6  (31) 

B7  (32) 

B8  (33) 
d = 4,430.30' 


W = 3. 8 kh  (6. 42  ft/sec) 
from  183°.  6 at  25'  elev 

M = 0.0058092 
e = 109°. 6 


T = 47.8°F  (8.8°C)  at  26'  elev 
C = 1, 104.8  ft/sec 


Mike 

True  Distance 
r2  (feet) 

True  Time 
r^c  (sec) 

Time  Data 
t(sec) 

Time  Error 
r^/c-t  (sec) 

1.) 

B6 

9,405.76 

8.514 

8.478 

.036 

2.) 

B7 

12,915.59 

11.690 

11.669 

.021 

3.) 

B8 

16,864.62 

15.265 

15.251 

.014 

True  (t^-t^)  = 3.177 
(tj't^)  Data  = 3.191 
True  = 6.751 

(tg-t^)  Data  = 6.773 
True  $2  = 3,509.83' 
TrueSg  = 7,458.86' 
Data  S.,  =3,525.42' 

sec 

sec  (long) 

sec 

sec  (long) 

Data  Sg  = 7,482.81' 
Wind  =3,516.79' 

1 “ 

( 

9,398.62' 

— 0.0759^  low 

2 

Wind  S =7,465.54'  ; 

1 wind 

9,556.96' 

1.614  high 
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[ 

Day:  12/3^4  (337) 

W6  (26) 

Event:  91 

W7  (34) 

Gun:  2 

W8  (37) 

t = 10:13:59.771  MST 
o 

d = 4,^30.33' 

t 

i- 

W = 3.8  kts  (6.42  ft/sec) 
from  183°.  6 at  25'  elev 

T = 47.8°F  (8.8°C)  at  26'  elev 
C = 1,104.8  ft/sec 

} 

f 

M = 0.0058092 

\ 

\ 

e = 199°. 6 

T. 

f ' 

i 

i 

True  Distance 

Mike  (feet) 

True  Time 
r^c  (sec) 

Time  Data 
t(sec) 

Time  Error 
(sec) 

[ 1.) 

W6  16,898.89 

15.296 

15.302 

^.006 

1 2.) 

W7  19,087.77 

\7.m 

17.269 

+.008 

i 3.) 

h 

W8  21,962.93 

19.880 

19.864 

+.016 

i! 

t 

t 

■ 

j 

'i 

true  (^2“*'])  ~ 1.981  sec 

(t2-t^)  Dato  = 1.967  sec  (short) 

True  (tg-tp  = 4.584  sec 

(tg-t^)  Data  = 4.562  sec  (short) 

True  $2  = 2,188.88' 

True  Sg  = 5,064.04' 

Trot 

%- 

^S-.'o 

i^sr!i 

= MS-fr 

1 

i 

Data  $2  =2,173.14'  1 

DataSg  = 5,040.10'  j 

= 16,789.74' 

0.646^  low 

1 

1 

1 

Wind  S = 2,148.89'  ) 

2 > R 

Wind  S'j  =4,991.61  ( wind 

= 16,988.98’ 

0.533^  KvaW 

Day: 

12/5/74  (339) 

B5  (30) 

d,  = 8,861. 

ir 

Event 

= 21 

B7  (32) 

2 

d- =4,430. 

36' 

Gun: 

2 

B8  (33) 

o 

D=  13,291. 

.47' 

" 12:48:00. 133  MST 

W = 

21.1  kts  (35.64  ft/sec) 

T = 56. 

,5°F  (13.6°C)at26 

' elev 

from  74**.  3 at  25'  elev 

C = 1/ 

114.2  ft/sec 

M = 

0.03198 

6 = 

218®. 9 

True  Distance 

True  Time 

Time  Data 

Time  Error 

Mike 

(feet) 

r^/c  (sec) 

t(sec) 

r/c-t  (sec) 

1.) 

B5 

7,026.36 

6.306 

6.204 

.102 

2.) 

B7 

12,915.59 

11.592 

11.438 

.154 

3.) 

B8 

16,864.62 

15.136 

14.880 

.256 

True  = 5. 286  sec 

^ (short) 

True  (tg-t^)  ~ 6.830  sec 
(ts-ti)  Data  = 8.676  sec  (short) 


True  - 

5,889.23' 

Tn.'eS^  = 

9,838.26' 

Data  $2  ~ 

5,831.72'  ] 

1 R = 8,947.46' 

Data  = 

9,666.80'  J 

I ^ 

Wind  = 
2 

5,611.18'  ] 

1 R = 10,307.47' 

Wind  Sg  = 

9,336.00'  ) 

wind 

Trut.  '^3=  7(1*1 


iliriiniMMIPliHRilRPlPI 


■ 1 \%  ir“  liirt 


Day:  12/5^4  (339) 
Event:  21 
Gun:  2 

t = 12:48:00. 133  MST 
o 


W = 21.1  kts  (35.64  ft/sec) 
from  74°. 3 at  25'  elev 

M = 0.03198 
0 = 128°.9 


187 


E3  (19) 

E2  (11) 

El  (8) 

d = 4,430.33' 


T = 56.5°F  (13.6°C)  at  26'  elev 
C = 1,114.2  ft/sec 


1.) 

2.) 

3.) 


Mike 

True  Distance 
(feet) 

True  Time 
r^/c  (sec) 

Time  Data 
t(sec) 

Time  Error 
ry'c-t  (sec) 

E3 

16,841.52 

15.115 

14.689 

.426 

E2 

19,023.42 

17.074 

16.607 

.467 

El 

21,895.15 

19.651 

19.068 

.583 

r 


True  (t^-t^)  = 1.958  sec 
(t2-ti)  Data  = 1.918  sec  (short) 
True  (tg-t^)  = 4.536  sec 
(tg-t^)  Data  = 4.379  sec  (short) 
True  $2  =2,181.90' 

True  Sg  = 5,053.63' 

Data  $2  = 2,137.04'  ) ^ 

Data  Sg  ~ 4,879.08'  | ® 

Wind  $2  “ 2,048.06' 

Wind  $2  = 4,701 


.06'  1 

.11'  i 


wind 


20,317.44' 

21,111.35' 


True  'Vj*  ‘hS-ro 

%=  n'? 

V,^,J  = ‘Mtl 


20.6^  high 
25.4'!^  high 
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Day:  12/5/74  (339) 
Event;  33 
Gun;  2 

t = 13:29:59.917  MST 


B6  (31) 

B7  (32) 

B8  (33) 
c<  = 4,430.30' 


1.) 

2.) 

3^ 


w = 

18.7  kts  (31.58  ft/sec) 
from  91°.  2 at  25'  elev 

T = 57.9°F  (14.4°C)at26' 
C = 1, 115.7  ft/sec 

elev 

M = 

0.02831 

e = 

202°. 0 

Mike 

True  Distance 
(feet) 

True  Time 
r^c  (sec) 

Time  Data 
t(sec) 

Time  Error 
r^c-t  (sec) 

B6 

9,405.76 

8.430 

8.386 

.044 

B7 

12,915.59 

11.576 

11.456 

.120 

B8 

16,864.62 

15.116 

14.874 

.242 

True  (t^-t^)  = 3. 146  sec 

(tj-t^ ) Date  = 3. 070  sec  (short) 

True  ~ 6.685  sec 

(*3“t|)  “ 6.488  sec  (short) 

True  $2  =3,509.83* 

True  $2  =7,458.86' 

Data  $2  =3,425.20'  ) ^ 

Data  $2  =7,238.66'  ) ® 

Wind  $2  =3,308.92'  ) ^ 

Winds  =7,006.11'  ) 


13,287.81' 


Trot. 


^ 41.3^  high 


4 


15,536.59' 


65.2^  high 


I I'J’Mhai'n 


-nfj^n‘rri^*^^T'*-^-'  T7-— -•• 
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Day:  ]2/5/74  (339) 
Event:  59 
Gun:  2 

t = 14:47:00.051  MST 
o 


A4(15)  82  = 8,860.63' 

A6(17)  82  = 4,430.59’ 

A7  (18)  0=  13,291.22’ 


\V  = 21.9  kts  (36.99  ft/sec)  T = 57.6°F  (14. 2°C)  at  26’ elev 
fr«n7»».5at25'ele»  c = M 15. 3 ff/».= 

M = 0.03316 
e = 146®. 3 


Mike 

True  Distance 
r2  (feet) 

True  Time 
r^c  (sec) 

Time  Data 
t(sec) 

Time  Error 
r^/c-t  (sec) 

A4 

6,983.02 

6.261 

5.668 

.593 

A6 

9,364.19 

8.396 

8.117 

.279 

A7 

12,885.60 

11.553 

11.260 

.293 

True  ~ ®®‘' 

(t2-t^)  Data  = Z449  sec  (long) 

True  (tg-t^)  = 5, 292  sec 

(t3-ti)  Data  = 5. 592  sec  (long) 

True  $2  = 2,381.17’ 

True  $2  = 5,902.58’’ 

Data  $2  = 2,731.37’  ) ^ 

Data  $2  = 6,236.76’  ) « 

Win8  =2,486.90’  ) 

2 > R 

Win8  $2  = 5,870.04  ) win8 


7,286.71’ 

7,877.07’ 


3\-S- 
= 33!? 


4.35^  high 


12.8^  high 
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Day:  12/5/74  (339) 
Event:  69 
Gun:  2 

14:47:00.051  MST 


N6  (6) 

N7(7) 

N8  (8) 

d = 4,430.34' 


w = 

21.9  kts  (36.99  ft/sec) 

T = 57. 

6°F  (14.2°C)at  26' 

elev 

from  79°. 5 at  25'  elev 

C = 1, 

115.3  ft/sec 

M = 

0.03316 

e = 

146°.3 

True  Distance 

True  Time 

Time  Data 

Time  Error 

Mike 

(feet) 

r^/c  (sec) 

t(sec) 

r^c-t  (sec) 

N6 

16,835.95 

15.095 

14.768 

.327 

N7 

19,020.91 

17.055 

16.697 

.358 

N8 

21,895.15 

19.632 

19.040 

.592 

True  (tj-t^)  = 1.959  sec 

(t^-t^)  Data  = 1.929  sec  (short) 

Tnje  (tg-t^)  = 4. 536  sec 

Data  = 4. 272  sec  (short) 

True  $2  =2,184.96' 

TrueSg  =5,059.20' 

Data  $2  “2,151.41'  ) ^ 

Data  Sg  =4,764.56'  ( o 

Wind  =2,029.17'  j 

2 >R' 

Winds.  =4,520.08'  ) wind 

w 


27,950.71' 

29,301.95' 


66. high 
74. 0^  high 


AO>A055  436  CALIFORNIA  UNIV  LOS  ANOELES  DEPT  OF  ATMOSPHERIC  SCIENCES  F/«  20/1 
PROBLEMS  AND  TECHNIQUES  OF  SOUND  RAN«IN6.(U) 

AUQ  77  M 6 WURTELE*  J ROE  DAAD07-75-C-0087 

UNCLASSIFIFD  Nt 

3aF3 


I 


END 


END 


8 -78  8-78 


Doy; 

12/5/74  (339) 

A5  (16) 

8^  = 4,430. 18' 

Event: 

83 

A6  (17) 

8^  = 8,860.79' 

Gun: 

2 

A8  (19) 

D=  13,290.97 

t = 

15:32:00.261  MST 

21.4  kfs  (36.  14  fj/sec) 
from  82°.  9 at  25'  elev 

0.03243 


T = 56.8  F (13.8  C)  at  26’  elev 
C = IJ  14.6  ft/sec 


True  Distance 


True  Time 


Time  Data 


Time  Error 


A5 

6,971.24 

6.254 

6.082 

.172 

A6 

9,364.19 

8.401 

8.093 

.308 

A8 

16,841.52 

15.110 

14.572 

.538 

True  (^2"*^])  “ 2. 147  sec 
(t2-t^)  Data  = 2.011  sec  (short) 
True  (tg-t^)  = 8. 855  sec 
(ts't^)  Data  = 8.490  sec  (short) 
True  $2  = 2,392.95' 

True  Sg  = 9,870.28' 

Data  $2  = 2,241.46'  I ^ 

Data  Sg  = 9,462.95'  ) ° 

Wind  “ 2, 1 17.43'  ) ^ 

Winds.  = 9,090.83'  I 


7,905.12' 

8,869.67' 


T'CUC.  ^3.1 


Day: 

12/5A4  (339) 

A6  (17) 

Event; 

: 85 

A7  (18) 

Gun: 

2 

A8  (19) 

t = 
o 

15:34:00.224  MST 

d = 4,430.39' 

W = 

22.8  kts  (38.51  ft/sec) 

T = 56.8°F  (13.8°C)at  26'  elev 

from  83°.4  at  25'  elev 

C = 1,114.6; 

Ft./sec 

M = 

0. 03455 

0 = 

150®. 2 

True  Distance 

True  Time 

Time  Data 

Time  Error  | 

Mike 

(feet) 

r^/c  (sec) 

t(sec) 

r^c-i  (sec) 

1.) 

A6 

9,364. 19 

8.401 

8.088 

.313 

2.) 

A7 

12,885.60 

11.561 

11.268 

.293 

3.) 

A8 

16,841.52 

15.110 

14.635 

.475 

True(f2-t^)  = 3. 159  sec 

Data  “ 3.180  sec  (long) 

True  = 6.709  sec 

^^3~^1^  Data  = 6.547  sec  (short) 

True  $2  = 3,521.41' 

True  $2  = 7,477.33' 

Data  $2  = 3,544.43'  ) ^ 

Data  Sg  = 7,297.29'  ) ® 

Wind  $2  = 3,411.61'  ) ^ 

Wind  S = 7,031.64'  ) wind 


26,704.98' 


30*1 
= 33.*3 


185^  high 
235^  high 


31,405.13' 
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r; 

f ; 


L* 


Day:  12/5/74  (339) 
Event:  85 
Gun:  2 

t = 15:34:00.224  MST 
o 


B5  (30) 

B6  (31) 

B7  (32) 
cl  = 4,430.55' 


I I ’ 


W = 22.8  kts  (38.51  ft/sec) 
from  83°.  4 at  25*  elev 

M = 0.03455 
6 = 209®.8 


T = 56.8°F  (13.8°C)  at  26'  elev 
C = 1,114.6  ft/sec 


Mike 

True  Distance 
(feet) 

True  Time 
r^c  (sec) 

Time  Dale 
t(sec) 

Time  Error 
r^/c-t  (sec) 

1.) 

B5 

7,026.36 

6.304 

6.356 

-.052 

2.) 

B6 

9,405.76 

8.439 

8.428 

+.011 

3.) 

B7 

12,915.59 

11.588 

11.431 

+.157 

True  ~ 2. 135  sec 

(t^-t^)  Data  = 2.072  sec  (short) 

True  (fg'fp  ~ 5.284  sec 

(t3-t^)  Data  = 5.075  sec  (short) 

True  $2  = 2,379.40' 

TrueSg  = 5,889.23' 

Data  $2  = 2,309.45'  1 ^ 

Data  $2  = 5,656.60'  ( ° 

Wind  $2  = 2,176.62  ) ^ 

Winds.  = 5,390.94'  I 


8,639.08' 

9,478.97' 


T'tuc  3\*.\ 
^3=  3 ft? 


22. H high 
34.9^  high 
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Day:  12/5/74  (339) 

Event:  87 

Gun:  2 

t = 15:36:00.237  MST 
o 

A4  (15) 

A6  (17) 

A7  (18) 

d2  = 8,860.63' 
d^  = 4,430.59’ 
D=  1 ^291.22* 

W = 

25.7  kts  (43.41  ft/sec) 

T = 56.8°F  (13. 8^0)  at  26'  elev 

from  77^.6  at  25'  elev 

C = 1 

,114.6  ft/sec 

M = 

0.03894 

0 = 

144°.  4 

True  Distance 

True  Time 

Time  Data 

Time  Error 

Mike 

Tj  (feet) 

r^c  (sec) 

t(sec) 

r^/c-t  (sec) 

A4 

6,983.02 

6.265 

5.626 

.639 

A6 

9,364.19 

8.401 

8.123 

.278 

A7 

12,885.60 

11.561 

11.225 

.336 

TrueCyt,)  = 2. 136  sec 

Dafa  = 2.497  sec  (long) 
TrueCr^-l,)  = 5.300  sec 

Data  = 5.599  sec  (long) 

True  $2  = 2,381.17’ 

True  $2  = 5,902.58' 

Data  Sj  ~ 2,783 
Data  Sg  ~ 6,240.65 
Wind  $2  = 2,502.59' 

Wind  $2  = 5,819.79* 


.16'  ) 
.65’  I 


wind 


7,638.95' 

8,329.02' 


%=  3o!<1 

3x:(. 


9.394  high 
19.34  higl> 


1 

Day:  12/5/74  (339) 

Event:  97 

Gun:  2 

t = 15:46:00.549  MST 
o 

E3  (19) 

E2(11) 

El  (8) 

d = 4,430.33' 

i 

W = 

22. 1 kts  (37.33  ft/sec) 

T = 56 

.7°F  (13.7°C)  at  26’  elev 

from  78°.  4 at  25’  elev 

C = 1, 

114.4  ft/sec 

M = 

0.03349 

! 

1 

0 = 

124°.  8 

1 

< 

True  Distance 

True  Time 

Time  Data 

Time  Error 

f 

Mike 

r2  {fe®0 

r^/c  (sec) 

t(sec) 

r^/c-t  (sec) 

1 

E3 

16,841.52 

15.113 

14.657 

.456 

2.) 

E2 

19,023.42 

17.071 

16.476 

.595 

3.) 

i 

1 

El 

21,895.15 

19.647 

18.884 

.763 

True  = 1.958  sec 

Data  = 1.819  sec  (short) 
True  (^3"^^)  “ 4.535  sec 
(t^-t^)  Data  = 4.227  sec  (short) 
True  $2  = 2,181.90' 

TrueSg  = 5,053.63' 

Data  Sj  ” 

Data  S. 


'H.i.i  = vfj 


Wind  $2  “ 
Wind  $2  “ 


= 2,027.09'  ) 
= 4,710.57'  ) 
= 1,942.40'  ) 
= 4,541.19'  j 


wind 


19,260.17' 

19,941.71' 


14.44  high 

18.44  high 
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Day:  } 2/5/74  (339) 

Event:  99 

Gun:  2 

t = 15:48:00. 215  MST 

0 

N5  (5) 

N7  (7) 

N8  (8) 

d2  = 8,860.67' 
dg  = 4,430.34' 
D=  13,291.01 

W = 

20.1  kts  (33.95  ft/sec) 

T = 56.7°F  (13.7®C)  at  26'  elev 

from  81  .Oat  25' elev  c=  1, 114.4  ft/sec 

M = 

0.03046 

e = 

147®. 8 

True  Distance 

True  Time 

Time  Data 

Time  Error 

Mike 

T2  (feet) 

T^c  (sec) 

t(sec) 

T^c-t  (sec) 

1.) 

N5 

15,632.00 

14.027 

13.705 

.322 

2.) 

N7 

19,020.91 

17.068 

16.570 

.498 

3.) 

N8 

21,895.15 

19.647 

19.141 

.506 

Tru«  (12“^])  “ 3.041  sec 
(t2“t^)  Data  = 2.865  sec  (short) 
True  (13“!^)  “ 5.620  sec 
(ta-ti)  Data  = 5.436  sec  (short) 
True  $2  = 3,388.91' 

True  $2  =6,263.15' 

Data  $2  =3,192.76'  ) 

DataS^  =6,057.88'  ) ° 

Wind  $2  =2,964 


Wind  $2  =5,715 


.35’  j 
.27’  ) 


wind 


14,769.15' 

15,527.35' 


Tv-v).  \fj=  If  S’*. 

4*,-  Vf.') 


^ 5.52i  low 


^ 0.67^  low 
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I 


Day:  \‘2/7/7A  (341) 
Event:  230 
Gun:  1 

t = 15:48:00.435  MST 


A3  (14) 

A2  (13) 

A1  (12) 

d = 4,430.70* 


W = 8. 4 kts  04. 19  ft/sec) 
from  354  .3  at  25*  elev 

M = 0.01277 
0 = 118°.9 


T = 53.4°F  (11.9°C)at26'  elev 
C = 1,110.9  ft/sec 


Mike 

True  Distance 
r^  (feet) 

True  Time 

Tj/c  (sec) 

Time  Data 
t(sec) 

Time  Error 
r^/c-f  (sec) 

1.) 

A3 

9,410.51 

8.471 

8.356 

.115 

2.) 

A2 

12,941.68 

11.650 

11.553 

.097 

3.) 

A1 

16,901.73 

15.214 

15.098 

.116 

True  “ 3. 179  sec 

(tj-t^)  Data  = 3. 197  sec  (long) 
True  (tg-t^)  = 6.743  sec 
(tg-t^)  Data  = 6.742  sec  (short) 
True  $2  =3,531.17* 

TrueS^  = 7,491.22* 

Data  = 3,551.55*  I ^ 

Data  Sg  = 7,489.69’  ( ® 

i Wind  $2  =3,524.20*  ) ^ 

Wind  =7,435.00*  I 


TVot  a 3^6 


10,856.16*  15.4qi:  high 

11,410.96*  — 21.31(  high 
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L. 


Doy:  1Z///74(341) 
Event:  230 
Gun:  1 


t = 15:48:00.435  MST 
o 


B6  (31) 

B7  (32) 

B8  (33) 
d = 4,430.30' 


W = 

8.4  kts  (14. 19  ft/sec) 

T = 

53.4°F  (11.9®C)  at  26* 

elev 

from  354®.  3 at  25*  elev 

C = 

1,110.9  ft/sec 

M = 

0.01277 

e = 

298®. 9 

True  Distance 

True  Time 

Time  Data 

Time  Error 

Mike 

Tj  (feet) 

r^/c  (sec' 

) t(sec) 

r^/c-t  (sec) 

1.) 

B6 

9,385.47 

8.449 

8.534 

-.085 

2.) 

B7 

12,891.24 

11.604 

11.700 

-.096 

3.) 

B8 

16,838.64 

15.158 

15.274 

-.116 

True  (^2"^])  ~ 3. 156  sec 


^^2”^1^  Data  = 3. 166  sec  (long) 


True  (tg-t^)  = 6. 709  sec 


(tg-t^)  Data  = 6.740  sec  (long) 


Tvrut  vV-j= 

^:i!C 

%wia  = ^\*  f 


True  $2  = 3,505.77* 


True  Sg  = 7,453. 17* 


Data  $2  = 3,517.11* 


Data  $2  = 7,487.47* 


Wind  $2  ~ 


Wind  S, 


= 3,544.45*  ) 
= 7,542.16*  j 


wind 


8,751.28* 

8,270.30* 


^ 6. 764  low 


11. 9<  low 
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9.  A MORE  GENERAL  VIEW  OF  THE  SOUND  RANGING  PROBLEM 

Die  foregoing  sections  show  censlusively  that  to  o greater  or  lesser  extent,  a 
number  of  commonly  occurring  circumstances  can  distxxt  the  acoustic  signal  of  the  Howitzer 
blast.  These  include  atmospheric  temperature  gradient,  on  absorbant  ground  surface,  and 
terrain  features  of  the  scole  of  hundreds  of  feet.  It  remains  true  that  in  theory  the  pulse 
front  of  the  signal  arrives  at  time  where  R is  the  shortest  distance  from  source  to 
microphone  and  c is  the  local  sound  speed.  However,  the  effect  of  diffraction  by  ony  of 
the  agencies  mentioned  above  is  to  deform  a relatively  sharp  pulse-front  --ideally  a step 
function — into  a signal  that  arrives  with  zero  amplitude  and  gradually  increases  to  the 
amplitude  given  by  the  theory  without  diffraction.  If  the  arriving  signal  hos  the  signature 
of  an  N-wave,  it  may  not  rise  to  this  amplitude  at  all. 

Thus  the  microphone  may  or  may  rK>t  be  cble  to  detect  the  true  arrival  time, 
depending  on  (1)  the  sensitivity  of  the  instrumentation,  and  (2)  the  background  noise 
level  at  the  time  of  arrival.  The  uncertainties  in  arrival  time  at  a number  of  microphones 
do  not  cancel  out,  but  may  lead  to  serious  rangir>g  errors. 

One  may  express  this  situation  in  a general  way  in  formulas.  If  the  signal 


is  an  undiffracted  pulse 


H(t  - H/c) 


P = P, 


and  is  detectable  at  any  two  similar  microphones.  Then 


- R|A  = ^2  " 


c(t2  - 1^) 


the  bosic  soond-ronging  assumption  (Section  I). 
However/  if  now 


P = P H(t  - H/c)  F(R,t) 
o 

where  F(R,t)  is  the  shape  factor  introduced  by  diffraction,  then  (9. 1)  no  longer  is 
valid.  We  now  equate  the  pressure  detectr^le  above  noise  level  by  the  tvro 
microphones,  and  obtain  an  equation  of  the  form 

F(R^,  t,  )=  F(R2,  tj).  (9,2) 

If  F(R,t)  were  expressible  as  a function  solely  of  the  single  variable  t - l^/c,  the  form 
of  F would  not  matter:  we  could  always  recover  (9. 1).  However,  the  signatures 
associated  with  diffroction  do  not  permit  this,  for  example,  that  given  by  (6.3). 

Thus  R|,  t^,  Rj,  and  t^  are  all  unknovm,  and  manipulations  and,  ultimately, 
approximations  are  required  to  obtain  the  desired  functional  form, 

R^-Rj  = ^(t^-tj)  . (9.3) 

To  produce  such  a formulation  might  not  be  impossible.  But  success  in  this  would 
imply  that  on  the  vrhole  the  diffraction  effects  were  smoll.  The  essential  feature 
of  a general  functional  equation  (9.2),  from  the  present  point  of  view,  is  that  the 
difference  tj  ” t^  does  not  uniquely  determine  the  difference  ^2  ~ basic 

sound-ranging  hypothesis  is  not  valid. 

This  situation  brings  us  face  to  foce  with  the  most  general  problem  of  sound 
ror^ir^.  A beginning  has  been  made  on  a general  approach  by  Miller  and  Engebos 
(1976),  but  it  is  not  clear  that  their  technique  can  accommodate  to  the  phenomenon 
of  strong  diffroction.  The  physics  of  the  problem  is  such  thot  observed  signol  orrivjl 


f 


time  differences  ore  not  in  generol  sufficient  information  to  determine  sourco 
location,  even  assuming  ideol  instrumentation.  If  assumptions  of  classical  sound 
parking  theory  are  not  valid  under  conditions  that  produce  diffroctive  effects, 
the  ranging  procedure  must  be  rethought  and  a much  more  complicated  informotion- 
processing  system  designed.  The  information  would  have  to  consist  of  the  time 
signatures  of  pressure  at  two  or  more  stations,  that  is,  amplitude  vs.  time  as  well 
as  time  of  arrival  information.  This  in  turn  means  coreful  calibration  of  the 
microphones.  Supposing  that  the  instrumentation  for  such  a system  existed.  How 
would  the  calculations  be  performed? 

First,  we  must  have  a furrctional  form  for  the  pressure  signatures,  derived 
from  a combination  of  theory,  computation,  and  observation.  Let  us  say 

p = F(R,  t - t^) 

where  t is  the  time  of  first  arrival  at  the  microphone,  unknown,  of  course.  Let 

Q 

microphone  number  k have  the  signature, 

pOc)  ^ ^ 

Then  for  each  microphone  there  are  two  unkrx>wns,  R and  t^.  Two  observations 
P.  at  times  t.  are  in  principle  adequate  to  determine  the  two  unkowns,  and  two 
microphones  are  sufficient  to  determine  the  source  locotion: 

p(k)  ^ F(F^\  t. -t^^  (9.4) 

I I Q 

In  practice,  as  is  usual  in  sound  ranging,  more  observations  will  be  useful:  the 
source  location  is  over-determined,  but  because  of  inaccuracies  in  practice,  a 


weighted  mean  of  these  estimates  is  likely  to  be  better  than  any  single  estimate. 
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The  soluHon  of  (9.4),  in  general  a transcendental  equation,  would  have 
to  be  done  by  a pre-programmed  calculator,  but  is  within  the  capabilities  of 
current  hand-held  equipment. 
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10.  CONCLUSIONS  AND  RECOMMENDATIONS 

(1)  The  exact  ranging  formulas,  without  approximation,  should  replace  the 
graphical  approximation  scheme  of  FM  6-122.  The  formulas,  including  correction 
for  constant  wind,  can  be  expressed  correct  to  the  first  order  in  the  Mach  number. 

These  formulas  ore  all  derived  from  first  principles  in  Section  1,  and  the  appropriate 
FORTRAN  softwear  is  provided  in  Annex  1. 

(2)  We  emphasize  that  these  formulas  referred  to  in  (1)  are,  or  soon  will  be,  within 
the  capability  of  hand-held  programmable  calculators,  and  recommend  that  this 
possibility  be  given  serious  consideration. 

(3)  The  Ho\/ltzer  blast  begins  as  a shock  front  and  rather  rapidly  decays  to  a sound 
pulse.  However,  both  theory  and  computation  (Section  7)  show  that  the  arrival 
time  of  the  pulse  at  sound-ranging  distances  is  measurably  in  advance  of  that  pre- 
dicted by  linear  theory.  This  suggests  the  form  of  a correction  term  in  the  sound 
rangirg  formulas.  In  Section  8 it  was  shown  that  the  incorporation  of  this  term 
improved  the  ranging  estimates  on  which  it  was  tested. 

It  is  recommended  that  further  tests  of  the  correction  term  (7. 12)  be  made 
on  independent  data  in  order  to  ascertain  whether  it  should  be  routinely  incorporated 
into  the  sound  ranging  formulas. 

(4)  In  Section  3 it  is  shown  that  the  character  of  the  soil  surface  significantly  lowers 
the  impedance  of  the  boundary  for  higher  frequencies  than  those  used  in  sound  ranging. 


ranging  errors. 

It  is  recommended  that  the  impedance  of  various  surfaces  be  determined 
for  sound-ranging  frequencies,  in  order  to  ascertain  whe  ther  this  effect  should 
be  irKorporated  into  the  procedure. 


(5)  The  theoretical  solution  of  a sound  pulse  propagating  in  on  atmosphere  with 
vertical  temperature  gradient  is  exhibited  in  Section  5.  Some  cases  were  numerically 
computed  and  it  was  concluded  that  diffraction  by  normal  atmospheric  temperature 
gradients  are  probably  not  a significant  source  of  error  in  sound  ranging. 

(6)  In  Section  6 it  is  shown  both  analytically  and  computationally  that  terrain 
features  of  the  scale  of  hundreds  of  meters  have  a profound  effect  on  the  form  of 
the  signal.  The  result  of  this  diffraction  is  that  the  signal  has  zero  amplitude  at  Hie 
time  of  arrival  and  increases  relatively  slowly  os  it  is  reinforced  continuously  by 
later  arrivals.  Depending  on  the  sensitivity  of  the  micropliones,  the  sound  ranging 
formulas  will  be  more  or  less  in  error. 

It  is  recommended  that  an  empirical  study  be  undertaken  to  determine  the 
threshhold  of  response  and  the  subsequent  record  of  the  microphone  exposed  to  a 
signal  with  the  characteristic  of  the  diffracted  wave.  With  Hiis  information,  together 
with  the  computations  of  this  report,  it  would  be  easy  to  verify  the  extent  of  the 
error  Introduced  by  terrain  features,  and  to  design  for  their  correction. 
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(7)  In  Section  8,  a subset  of  the  PASS  data  was  examined.  The  atmospheric 
conditions  in  those  events  selected  were  most  nearly  conformir^g  to  the  assumptions 
underlying  sound  ranging  technique.  Although  the  ranging  formulas  were  highly 
successful  in  target  acquisition  in  certain  of  these  cases,  in  some  others  the  ranging 
errors  were  off  by  as  much  os  50  to  100  percent.  Using  different  sets  of  microphones, 
we  obtained  both  excellent  and  poor  ranging  estimates  in  the  same  event.  After 
eliminating  other  sources  of  error,  we  concluded  that  the  source  was  instrumental, 
most  likely  in  the  timing  system. 

It  is  recommended  that  a controlled  experiment  be  designed  to  test  the 
quantitative  reliability  of  the  instrumentation  in  the  field,  in  order  to  obtain  a 
better  understanding  of  the  source  of  ranging  errors. 

(8)  This  report  has  not  been  able  to  treat  the  problem  of  diffraction  by  wind  shear. 

It  is  clear,  however,  that  this  effect  con  be  of  the  same  importance  as  diffraction 
by  terrain  features,  but  that  the  theoretical  methods  of  dealing  with  it  ore  quite 
different.  It  is  therefore  recommended  that  diffraction  by  wlndsheor  be  investigated 
in  relation  to  the  problem  of  sound  ranging. 

(9)  In  Section  9 it  is  pointed  out  that  the  use  of  further  information  (that  Is,  other 
than  arrival  time)  from  the  microphone  records  could  improve  ranging  estimates. 

This  point  is  also  made  by  W.S.M.R.  personnel  Miller  and  Engebos  (1976). 

It  is  highly  recommended  that  the  possibility  of  obtaining  and  using  quantitative 
signal  information  be  investigated.  This  research  effort  would  necessaril.  Involve 


the  interaction  of  expertise  in  instrumentation,  in  atmospheric  physics,  and  in 
mathematics;  but  the  chief  hope  of  improvement  in  sound  ranging  lies  in  this 


direction. 


11.  REFERENCES 


Brekhovskikh,  L.  M. , 1960:  Waves  in  loyered  Medio,  Academic  Press. 

Brode,  H.,  1959;  "Blast  wave  from  a spherical  charge,"  Phys.  Fluids,  2,  217-229. 

Campbell,  R.  G. , 1958:  "Initial  wave  phenomena  in  a weak  spherical  blast," 

J.  Appl.  Phys. , 29,  55-60. 

Carlson,  M. , and  Hunt,  J.  T. , 1974:  "Effect  of  geometric  symmetry  on  acoustic 
scattering  in  the  shadow  region,"  J.  Acoust.  Soc.  Am. , 55,  964-967. 

Chessell,  C.  I. , 1976:  "Observations  of  the  effects  of  atmospheric  turbulence  on 
low-frequency  sound  propagation,"  J.  Acoust.  Soc.  Am. , 60,  29-33. 

Chien,  C.  F. , and  Soroka,  W.  W. , 1975:  "Sound  propagation  along  an 
impedance  plane, " J.  Sound  Vib. , 43,  9-20. 

Delaney,  M.  E. , and  Bazley,  E.  N. , 1970:  "Monopole  radiation  in  the  presence 
of  an  absorbing  plane,"  J.  Sound  Vib. , 13  , 269-279. 

N.  , , 

Dickenson,  P.  J, , and  I>oak,  P.  E. , 1970:  "Measurements  of  the  normal  acoustic 
impedance  of  ground  surfaces,"  J.  Sound  Vib.,  13  , 309-322. 

Dock,  P.  E. , 1952:  "The  reflexion  of  a spherical  acoustical  pulse  by  an  absorbant 
plane,"  Proc.  Royal  Soc. , A,  215,  233-254. 

Donato,  R.  J. , 1960b:  "Spherical  vrave  reflection  from  a boundary  of  reactive 

impedance  using  a modification  of  Cagniard's  method, " J.  Acost.  Soc.  Am. 

Donato,  R.  J. , 1976a:  "Propagation  of  a spherical  wave  near  a plane  boundary 
with  complex  impedance, " J.  Acost.  Soc.  Am. , 60,  34-39. 

Bnbleton,  T.  F.  W. , 1956:  "The  propagation  and  reflection  of  sound  pulses  of 
finite  amplitude,"  Proc.  Phys.  Soc.  London,  B69,  382-395. 

Embleton,  T.  F.  W. , and  Piercy,  J.  E. , 1974;  "Effect  of  ground  on  near-horizontal 
sound  propagation,"  Society  of  Automotive  Engineers,  Congress,  Detroit, 
Michigan,  1974. 


i 


Embleton,  T.  F.  W. , Piercy,  J.  T. , and  Olson,  N. , 1976:  "Outdoor  sound 

propagation  over  ground  of  finite  impedance,"  J.  Acoust.  Soc.  Am.,  59,  267-277, 


208 


I 


Embleton,  T.  F.  W. , Thiessen,  G.  T. , and  Piercy,  J.  E. , 1976:  "Propagation 
in  on  inversion  and  reflections  at  the  ground, " J.  Acost.  Soc.  Am. , 59, 

278-282. 

Field  Manual  FM  6-122,  1964:  "Artillery  Sound  Ranging  and  Flash  Ranging", 

Department  of  the  Army,  Washington,  D.C. 

Franck,  H. , and  Sager,  G. , 1975:  "Zur  Schallv/egverteilung  in  der  unteren 
Grundschicht  bei  Vorhandsein  einer  Bodeninversion,  " Z.  fur  Meteor. , 

^ 130-137. 

Friedlander,  F.  G.,  1958:  Sound  Pulses,  Cambridge  University  Press. 

Huang,  H. , 1975:  "Scattering  of  spherical  pulses  by  a hard  cylinder," 

J.  Acoust.  Soc.  Am. , 58,  310-317. 

Ingard,  U. , 1951:  "On  the  reflection  of  a spherical  sound  wave  from  an  infinite  plane," 
J.  Acost.  Soc.  Am.,  23  , 329-335. 

Kinney,  G.  F. , 1962:  Explosive  Shocks  in  Air,  Macmillan  Co. 

Krasilnikov,  V.  A.,  1963:  Sound  and  ultrosound  waves  in  air,  water,  and  solid  bodies. 
Published  for  the  National  Science  Foundation  by  the  Israel  Program  for 
Scientific  Translation,  Jerusalem. 

Lov^ead,  R.  B. , and  Rudnick,  I. , 1951:  "Measurements  of  an  acoustic  wave 
propagated  along  a boundary, " J.  Acost.  Soc.  Am. , 23,  541-549. 

Lee,  E.  L. , Hoinig,  H.  C. , and  Kury,  J.  W. , 1968:  "Adiabatic  expansion  of 
high  explosive  detonation  products,"  UCRL-50422,  Lawrence  Radiation 
Laboratory,  Univ.  of  California. 

McNicholos,  J.  V.,  Uberall,  H. , and  Choate,  K.,  1968:  "Pulse  shapes  of 

creepirig  waves  around  soft  cylinders,”  J,  Acost.  Soc.  Am. , 44,  752-764. 

Meecham,  W.  C. , 1956:  "On  the  use  of  the  Kirchhoff  approximation  for  the 

solution  of  reflection  problems",  J.  Rotionol  Mech.  Anol. , 5,  323-334. 

Meecham,  W.  C. , 1957:  "Vorational  method  for  the  calculation  of  the  distribution 

of  energy  reflected  from  a periodic  surface,"  J.  Acost.  Soc,  Am.,  27,  361-367. 

Miles,  J.  M. , 1967:  "Decay  of  spherical  blost  waves, " Phys.  Fluids,  10,  2707-2708. 


I 


J 


209 


Miller,  W. , and  Engebos,  B. , 1976;  “A  mathemafical  sfruch/re  for  refinement 
of  sound  ranging  estimates,”  R and  D Tech  Rept.  ECOM-5805,  ASL, 

White  Sands  Missile  Range. 

Milne,  E.  A.,  1921:  "Sound  waves  in  the  atmosphere,"  Phil.  Mag. , 6th  ser. , 

42,  96-114. 

Morse,  P.  M. , and  Ingard,  U. , 1968:  Theoretical  Acoustics,  McGraw-Hill 
Book  Co. 

Naugolnykh,  K.  A. , 1966;  "Compression  wave  radiated  by  an  expanding  sphere, " 
Soviet  Physics-Acoustics,  11,  296-301. 

Piercey,  J.  E. , and  Embleton,  T.  F.  W. , 1974:  " Effect  of  weather  and  topography 

on  the  propagation  of  noise  at  Vancouver  Airport, " National  Research 
Council  of  Canada,  Division  of  Physics,  Acoustics  Section. 

Proud,  J.  M. , Tamarkin,  P.,  and  Meecham,  W.  C.,  1958:  "Reflection  of  sound 
from  a surface  of  saw-tooth  profile,"  J.  App.  Phys. , 28,  1298-1301. 

Rudgers,  A.  J. , and  Uberall,  H. , 1970:  "Pulses  specularly  reflected  by  a soft 
cylinder,"  J.  Acoust.  Soc.  Am. , 48,  3710381. 

Rudnick,  I. , 1947;  "The  propagation  of  an  acoustic  wove  along  a boundary, " 

J.  Acost.  Soc.  Am.,  19,  348-356. 

Sedov,  L.  I.,  1959;  Similarity  and  dimensionol  methods  in  mechanics.  Academic 
Press. 

Taylor,  G.  I,,  1946;  "The  air  wave  surrounding  an  expanding  sphere," 

Proc.  Royal  Soc.,  A,  186,  273-292. 

Taylor,  G.  I. , 1950;  "The  dynamics  of  the  combustion  products  behind  plone 

and  spherical  detonation  fronts  in  explosives, " tVoc.  Royal  Soc. , A,  200, 
235-247. 

Taylor,  G.  I. , 1950:  "The  formation  of  a blast  wave  by  a very  intense  explosion, " 
Proc.  Roy.  Soc.  London,  A 201,  159-186. 

Trulio,  J.  G.,  Perl,  N.  K.,  and  Balanis,  G.  N.,  1976;  "The  Afton  2A  computer 
code  revised  user's  manual,  AFWL-TR-75-111,  Defense  Nuclear  Agency, 
A.F.  Weopons  Loborotory,  Klrtland  A.F.B. , N.M. 


210 


i 


I ‘ 
i 

I 


I 

( 


I 

I 


Uberall,  H. , Doolittle,  R.  D. , and  Me  Nicholas,  J.  V.,  1966:  "Use  of  sound 
pulses  for  a study  of  circumferential  waves,"  J.  Acoust.  Soc.  Am. , 39, 
564-578. 

Thomosson,  S.  -I. , 1976:  "Reflection  of  waves  from  a point  source  by  an  impedance 
boundary",  J.  Acost.  Soc.  Am. , 59,  780-785. 

U.  S.  Army  Field  Artillery  School,  1972:  "Development  of  the  Field  Artillery 
Observation  (Target  Acquisition)  Battalions,"  Fort  Sill,  Oklahoma. 

Wenzel,  Alan  R. , 1974:  "Propagation  of  waves  along  an  impedarKe  boundary", 

J.  Acost.  Soc.  Am. , 55,  956-963. 

Whithorn,  G.  B. , 1974:  Linear  ond  Nonlinear  Waves,  John  Wiley  and  Sons. 

Whithorn,  G.  B. , 1950:  "The  propagation  of  spherical  blast,"  Proc.  Royal  Soc. 

A,  203,  571-581. 

Wurtele,  M.  G. , 1961:  "On  the  problem  of  truncation  error,"  Tellus,  13,  379-391. 


! 

I 


i 


I 

I 

1 

i 

i 

1 

I I 


J 


i 


ANNEX  I 

FORTRAN  CODE  FOR  SOUND  RANGING  ESTIMATES 
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